Volume 7 


APRIL, 1936 


PHYSICS 





A Journal of General and Applied Physics 





The Dielectric Constant of Air at Radiofrequencies 


L. G. Hector AnD H. L. Scnuttz, Department of Physics, University of Buffalo 
(Received November 16, 1935) 


The heterodyne beat frequency method of determining 
dielectric constants has been applied to the measurement of 
the dielectric constant of air with an applied frequency of 
field of 900,000 cycles per second. Radio broadcast signals, 
stripped of modulation, were used as sources of alternating 
current of constant high frequency. An electron coupled 
oscillator was used in the circuit containing the test 
condenser. Both sources of high frequency were isolated 


from the detector through buffer stages of amplification. 
Determination of the beat frequencies was made by 
photographic recording with a two string oscillograph. A 
dielectric constant of 1.00058986 +0.00000050 is found for 
air at normal pressure and temperature. This result is in 
close agreement with the results of other investigators for 
lower frequencies of applied field and for static measure- 
ments. 





INTRODUCTION 


N the course of preliminary work on the effects 
of intense electric fields on the dielectric con- 
stants of gases, we have applied the heterodyne 
beat frequency method to the measurement of 
the dielectric constant of air when the frequency 
of alternation of the applied field is 900,000 
cycles per second. 

Basically this method of measuring dielectric 
constants is old. It consists in mixing a current 
of high but unvarying frequency in a detector 
with current supplied from a variable frequency 
oscillator. The detector then delivers a current 
with the difference or beat frequency to some 
type of frequency indicating instrument. 

The variable frequency oscillator circuit in- 
cludes a condenser whose capacity may be 
changed by the withdrawal or insertion of the 
material whose dielectric constant is to be de- 
termined. This capacity change may be com- 
pensated for by altering the setting of another 
condenser in the circuit; or the changes in 
capacity may be computed from the change in 
- frequency of the beat note. 

The principal experimental difficulties with the 
method are : 


. Maintaining the desired constancy of frequency of a 
standard high frequency source. 

. Maintaining the constancy of a variable frequency 
oscillator with respect to causes other than the change in 
the capacity of the test condenser. 

3. Avoiding the tendency of the two sources to synchronize 
when the beat frequency approaches zero. 


A solution to the first of these difficulties is the 
use of a crystal controlled oscillator with tem- 
perature control for the crystal and this type of 
source for a standard high frequency current is 
now in general use by experimenters in this field. 

The third difficulty, that of synchronization, 
may be handled by radiation coupling to the 
detector from the two sources of high frequency 
current as was reported by Bryan and Sanders! 
in 1928. 

It can also be solved by the use of buffer 
stages of.amplification such as have now come to 
be accepted technique in high frequency radio 
engineering. 

Solution to the second difficulty, that of the 
constancy of the variable frequency oscillator to 
causes other than changes in the test condenser, 
is now possible to a greater extent than formerly 


1A. B. Bryan and I. C. Sanders, Phys. Rev. 32, 302 
(1928). 
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through the development of multigrid tubes, 
special oscillating circuits for use with such 
tubes, and buffer stages to isolate the oscillator 
from the effects of any possible load variations 
in the output circuit. 


EXPERIMENTAL METHOD 


We have attempted to incorporate the im- 
provements outlined above in this well-recog- 
nized method for the measurement of dielectric 
constants with the hope that sufficiently im- 
proved performance might be obtained to enable 
us to look for dielectric saturation phenomena in 
gases. Largely as an attempt to see what pre- 
cision could be attained we have made a series 
of measurements on the dielectric constant of 
air and wish to report the results of these meas- 
urements as well as the details of the method in 
the present paper. 

Fig. 1.shows the circuit diagram of the variable 
frequency oscillator and its buffer stage. The 
oscillator is the so-called electron coupled type. 
Essentially it is a Hartley oscillator where the 
grid of the tube used ordinarily as a screen is 
used as the plate. Output power, however, is 
taken from the actual plate of the tube. An 
analysis of this type of oscillator and the reasons 
to which its constancy of frequency may be 
attributed have been given at length in papers 
by F. B. Llewellyn’ and by J. B. Dow.’ 

The oscillator is resistance-capacity coupled to 
a buffer stage which in turn supplies power to 
the detector circuit. The buffer stage has its 
grid biased for class B operation and the output 


2? F. B. Llewellyn, Proc. I. R. E. 19, 2063 (1931). 
3]. B. Dow, Proc. I. R. E. 19, 2095 (1931). 
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circuit may be tuned either to the fundamental 
or to an harmonic of the oscillator. In the work 
reported in this paper the second harmonic 
was used. 

The entire oscillator unit was built into a 
shield box made of 3/8 inch sheet aluminum. 
Such heavy construction was found necessary to 
insure the required mechanical and hence elec- 
trical capacitative stability. The buffer stage 
was built into a separate aluminum shield and 
the two were kept in a large wood box to improve 
temperature control. 

Filament and plate batteries were shielded in 
a copper lined box. 

In Fig. 1 the test condenser is C. C; is a vernier 
condenser used to set the oscillator at approxi- 
mately the frequency desired. C, is a smaller 
vernier condenser used for fine adjustment of 
the oscillator to zero beat with the standard 
high frequency. 

The test condenser was of the multiple plate 
parallel plate type. The plates were made of 
brass and the whole unit was housed in a con- 
tainer turned from solid brass. 

The source of standard high frequency was the 
broadcast signal from a radio station (WBEN) 
equipped with late type temperature controlled 
crystal oscillators. A method for removing the 
modulation from such signals so that the presence 
of program need not interrupt the use of the 
carrier for laboratory measurements is described 
in another paper by the authors. 

The two high frequencies were mixed in the 
detector amplifier circuit shown in Fig. 2. The 
output of this circuit feeds a standard string 
oscillograph. 

Some idea of the constancy both of this 
standard signal and of the local oscillator may 
be gleaned from frequency data of which the 
following is typical. The beat note between the 
second harmonics of the standard frequency 
(900,000) and that of the local oscillator was 
set at about 100 cycles and the actual beat 
frequency accurately determined at 20-minute 
intervals. The values found were 99.94, 99.69, 
99.00 showing a variation of about 1 cycle per 
second in 1,800,000 cycles per second in 40 
minutes. 


4L. G. Hector and H. L. Schultz, Rev. Sci. Inst. 7, 139 
(1936). 
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EXPERIMENTAL PROCEDURE 


Before taking any data the oscillator was 
allowed to run for a period of ten to twelve 
hours. The test condenser was then evacuated 
and the local oscillator set so that its second 
harmonic (1,800,000) gave zero beat with that 
of the second harmonic of the standard signal. 
This condition could be observed by the motion 
of the needle of a d.c. milliammeter in the plate 
circuit of the amplifier connected to the de- 
tector (see Fig. 2). 

Air was now allowed to enter the test con- 
denser after passing through a solution of KOH 
for removing COs, and through two large reser- 
voirs packed with CaCl, for removing water 
vapor. For the change from vacuum to atmos- 
pheric pressure, the change in the frequency of 
the oscillator was approximately 215 cycles per 
second. This corresponds to a frequency of about 
430 cycles per second between the second har- 
monic frequencies. 

The actual frequency was determined pre- 
cisely by means of photographic recording with 
a two string oscillograph where the beat note was 
recorded by one galvanometer and a timing 
wave by the other. The timing employed was 
that of the 60-cycle electric power supply. 

Checks of the frequency deviations of this 
source from 60 cycles per second show that the 
variations seldom exceed 0.05 cycle per second 
and that they are as often above as below the 
assigned value. 


DETERMINATION OF STRAY AND EFFECTIVE 
CAPACITIES IN THE CIRCUIT 


In computing the change in capacity of the 
test condenser from beat frequency data it 
becomes necessary to know the ratio of the 
stray capacity in the oscillatory circuit to that 


FROM 
BUFFER 


59 OPERATED AS TRIODE 
TO OSCILLO- 
GRAPH OR 
SPEAKER ‘ 


*280V 








Fic. 2. Detector and A.F. amplifier. 


CONSTANT 


OF AIR 135 
of the effective capacity of the test condenser. 
This necessity arises from the fact that ob- 
viously the stray capacity of the circuit does 
not change when gas is introduced into the test 
condenser. 

The capacity of the test condenser was first 
determined by means of substituting a new, 
calibrated, type 222F General Radio air con- 
denser in the circuit. The mean of several 
determinations gave a value of 518.4 mmfd. 
It was estimated that as much as 1.5 mmfd. 
of this capacity was so located as not to be 
effected by the introduction of gas. Hence the 
effective value of the test capacity is 516.9 mmfd. 

For determining the ratio of stray to effective 
capacity for the entire circuit the following 
procedure was employed. 

As in the experimental procedure described 
above, the frequency change of the variable 
frequency oscillator, originally at 900,000 cycles 
per second, was observed for an air pressure 
change in the test condenser from vacuum to a 
value near atmospheric. 

Capacity from a variable standard condenser 
was now added in parallel to the test condenser 
until the frequency with the test condenser 
evacuated was 700,000 cycles per second. (Stand- 
ard frequency was obtained from the signal of 
radio station WLW.) The frequency change on 
filling the test condenser with air to the same 
pressure as before was then determined. 

From a knowledge of the two original oscillator 
frequencies, the two beat frequencies, the 
effective value of the test condenser, and the 
amount of added capacity for the second series 
of data the desired ratio between the stray 
capacity of the circuit and effective capacity 
of the test condenser, R, can be determined. The 
value for our circuit was found to be 0.09189. 


CALCULATION OF THE DIELECTRIC CONSTANT 


If we write ¢e for the dielectric constant of the 
gas under the conditions of pressure and tem- 
perature at which the measurement is made, 
we can compute the value of (e—1) from a 
knowledge of the ratio, R, and the change in 
frequency of the oscillating circuit when the 
surroundings of the test condenser are changed 
from vacuum to this pressure of gas. 
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Equations showing this relation may be de- 
veloped as follows. The natural frequency, f, of 
the oscillating circuit is given by 


f=1/2xr(L(C+c))}, (1) 


where L is the inductance, C the effective value 
of the test condenser, and c the stray capacity 
of the circuit. 


df /dC = —1/44L\(C+0)!. 


Remembering that dC=C(e—1) and_ that 


R=c/C, we may write 
e—1=—2df(1+R)/f. (2) 


Since the change in frequency is rather large 
to be treated in the manner of a differential, a 
more precise expression for (e—1) may be ob- 
tained by expanding (1) by the method of 
increments. The result obtained in this manner is 


e—1=(14+R)[1/(1+4f/f2—-1]. (3) 


Results published in this paper have been 
computed with the aid of Eq. (3) although the 
differences as compared with those obtained with 
Eq. (2) are not large. 


RESULTS 


Table I shows the results of 26 determinations 
of the dielectric constant of air corrected to 
normal temperature and pressure. 

It shows the smallness of errors that are 
accidental in character. However, constant 
sources of error decrease the reliability of the 
result. The most important of these lies in the 
determination of the ratio of the stray to the 
effective capacity in the circuit. An error of 
1 percent in this ratio will affect the final results 
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TABLE |. Experimental values for dielectric constant of air 
at N. ar 








1.0005895 
1.0005892 
1.0005900 
1.0005895 
1.0005897 
1.0005902 
1.0005901 
1.0005900 
1.0005899 


1.00058986 
Average deviation 


from the mean 0.00000043 
Number of observations 26 


1.0005899 
1.0005899 
1.0005897 
1.0005893 
1.0005902 
1.0005903 
1.0005903 
1.0005907 


1.0005899 
1.0005895 
1.0005902 
1.0005910 
1.0005888 
1.0005903 
1.0005888 
1.0005889 
1.0005906 
Average 








for (e—1) by 0.08 percent and this error is 
probably the limiting factor in the reliability of 
our data. 

It may be noted that neither Eq. (2) or (3) 
contains a term for either the stray or the 
effective capacity of the circuit independently 
but only one for the ratio of these two quantities. 
Precision of measurement of these quantities is 
then only important insofar as it effects this 
ratio and an error in this ratio affects the value 
of (e—1) only to the extent indicated above. 

On this basis the value for the dielectric 
constant of air at normal pressure and 0°C 
(when the air is dried with CaCl. only and has 
had CO. removed by a solution of KOH) at 
900,000 cycles per second is 


1.00058986 
+0.00000050. 


This value compares with 1.000591 found by 
Jordan, Broxon, and Walz® at 70,000 cycles per 
second in 1934 and with 1.000590 found by 
Boltzmann by static methods as long ago as 1874. 


>A. R. Jordan, Jamies W. Broxon and Frank C. Walz, 
Phys. Rev. 46, 66 (1934). 





APRIL, 1936 


PHYSICS 


VOLUME 7 


Low Pressure Arc Characteristics* 


J. D. Cosine, Graduate School of Engineering, Harvard University 
(Received July 29, 1935) 


The results of an oscillographic study of the short a.c. arc 
between pure graphite electrodes show that the reignition 
potential is a linear function of the gas pressure for air, 
No, Oz, and COs. The relation is expressed by V-=A+Bp. 
The constant A follows the minimum sparking potential of 
the gas, increases with the rate of rise of the impressed 
voltage, and is independent of the arc current, except for 


CO,. The constant B depends on the gas, arc current, and 
rate of rise of impressed voltage. The reignition potential 
is closely related to the glow discharge and is shown to 
approach the sparking potential of the gas as the r.m.s. 
value of the current approaches zero. The factors that affect 
reignition are indicated, and a reasonable mechanism of the 
phenomenon is stated. 





INTRODUCTION 


N an experimental study of the factors affect- 

ing the reignition of a short a.c. arc between 
refractory electrodes, it has been shown! that for 
air the reignition potential increases with pres- 
sure, decreases with increasing current, and is 
markedly reduced by the addition of small 
amounts of low work-function impurities. The 
experimental results presented in this report in- 
dicate that for several gases there is a linear 
relationship between the reignition potential and 
the gas pressure in the case of a short arc between 
refractory electrodes. 


EXPERIMENTAL RESULTS 


The apparatus and circuit are the same as 
previously described.! The power was supplied to 
the arc by a power transformer, and the r.m.s. 
arc current was limited to the desired value by a 
resistance in series with the arc and the trans- 
former. Oscillograms were taken of the voltage 
across the gap by means of a magnetic oscillo- 
graph. The discharge investigated took place 
between the plane end surfaces of cylindrical 
electrodes of 6.4 mm diameter separated by a 
gap of 1 mm. The electrodes were placed under a 
large bell jar which was connected to suitable 
pumps and gauges. Data were taken for pressures 
above 10 cm of Hg. The results for pressures 
below this were not believed to be representative 
of the electrodes being studied, for the discharge 
was not confined to the electrodes; but covered 
all exposed parts of the electrode holders. 

* Part of the experimental results of a thesis in partial 
fulfillment of the requirements for the degree of Doctor of 
Philosophy, California Institute of Technology, June 1934. 


1S. S. Mackeown, J. D. Cobine and F. W. Bowden, 
Elect. Eng. 53, 1081 (1934). 


The electrodes were spectrographic carbons 
supplied by the National Carbon Company. Ac- 
cording to the manufacturer these carbons were 
made of specially purified graphite having con- 
siderably less than 0.001 percent ash; examina- 
tion had show the presence of silicon, calcium, 
magnesium, silver, and copper. It was stated that 
no carbon examined had shown the presence of 
more than two or three of these elements and 
some had shown the lines of carbon alone. As was 
pointed out previously,' electrodes of high purity 
exhibited very uniform reignition potentials from 
cycle to cycle. 

The gases used were nitrogen, oxygen, and 
carbon dioxide of commercial purity, and air. 
The electrodes were cleaned and the gap adjusted 
at frequent intervals so that the degree of 
uniformity of surface and gap would be as high 
as possible. Each value of reignition potential is 
the average of about the first forty individual 
values taken from oscillograms. Deviations from 
the average were small. In the case of nitrogen it 
was found that no error was introduced when the 
deflections were taken directly from the oscillo- 
graph screen, that is, the time studied could 
cover several seconds. 

There was practically no destruction of the 
electrodes when nitrogen was used. With oxygen 
there was considerable consumption of the 
graphite; an after-glow of short duration was 
noted between the electrodes when a current 
greater than 4 amperes was interrupted, indicat- 
ing, possibly, sustained combustion. This effect 
was not noted in air, neither was the electrode 
consumption very large. The use of carbon 
dioxide resulted in considerable burning of the 
electrodes, probably due to the fact that COsz is 
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Fic. 1. Arc reignition potential for 
graphite in nitrogen. Transformer 
voltage=2200. Parameter=r.m.s. 
are current. 


voltage = 2200. 


are current. 


dissociated at high temperatures into O and CO. 

The are current is given in r.m.s. amperes, as 
indicated by an iron vane a.c. meter. Needless to 
say, the current wave is not sinusoidal; but has a 
very low or zero value during the period from arc 
extinction to reignition. 

The experimental results showing the reigni- 
tion potential as a function of gas pressure for 
different currents and gases are presented by the 
curves of Figs. 1 to 3. It is seen that in each 
case the reignition potential of the arc is a linear 
function of the pressure over a wide range of 
pressure; following a curve such as V,=A+ Bp. 
The dotted extension of the curves indicate the 
value of A. It should be noted that, except for 


TABLE I. Reignition constants for graphite. 








TRANS- 
FORMER A 
VOLTAGE (volts) 


2200 275 
3300 325 


MINIMUM Bi 
SPARKING (volts/cm Hg) 
POTENTIAL? lamp. 6.5 amp. 


251-298 58.7 35.0 
27.0 


GAs 


Nitrogen 





Air 2200 
3300 


325 327 
375 


40.0 
25.0 


2200 
3300 


Carbon 
Dioxide 


420 25.0 


425 to 600° 
. 25.0 


400 to 600 
2200 550 


Oxygen 439-455 


25.0 











! Determined for the straight line portion of the curves. 

? Range for currents 11 to 1 amperes. 

*J. J. Thomson, Conduction of Electricity in Gases, Vol. 2, third edi- 
tion, p. 487. 


Fic. 2. Arc reignition potential for 
graphite in oxygen. Transformer 
Parameter=r.m.s. 








0 3 
Pressure incm of Hg 


Fic. 3. Arc reignition potential for 
graphite in carbon dioxide. Trans- 
former voltage=2200. Parameter 
=r.m.s. arc current. 


carbon dioxide, each curve of the family has the 
same intercept at p=0 for a given gas. The linear 
relation has been found to extend to atmospheric 
pressures for Oz, Ne, and air. In this case a higher 
transformer voltage was required so that the 
results may not be compared directly. Table I 
summarizes the constants for the two transformer 
voltages indicated. 

The curves of Figs. 4 to 6 show the dependence 
of V,, at constant gas pressure, upon the arc 
current. With nitrogen and with carbon dioxide 
there are two nonlinear regions in the range of 
current covered, just as was shown! to be true for 
air. One region for low currents, up to a little 
more than 6 amperes for N»2 and up to about 4 
amperes for CO.; and the other region involving 
higher currents. However in the case of oxygen 
the reignition potential decreases almost linearly 
with increasing current. The curves of Fig. 4 
indicate that nitrogen as well as air! has, under 
the conditions of these experiments, a current at 
which the curves of constant pressure cross. 
When the current is increased above this value 
the highest pressures have the lowest V,. This 
critical value of current marks the limit of the 
linear relationship between V, and pressure for 
air and nitrogen. This point was not reached for 
Oz and COs. In all these curves the sparking 
potential is indicated for each pressure at zero 
current. The V,—J curves have the same shape 
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Fic. 4. Arc reignition potential for 
graphite in nitrogen as a function of 
the arc current. Parameter is gas 
pressure in centimeters of mercury. 


at the higher transformer voltage as for those 
shown in the figures. 

Fig. 7 shows for oxygen and for nitrogen the 
dependence of the constant B upon the arc 


current. For oxygen this constant is nearly a 
linear function of the arc current. For nitrogen 
the constant B was found to depend upon J, J°, 
and J* as shown by the equation: B=67.7 —9.9] 
+1.35/?—0.1J*. For nitrogen the critical value of 


B -Vo/ts per Centimeter of Mercury 


! 4 678 9 1 
re Current inAmperes 


Fic. 7. The constant B as a function of the are current. 
A—circles, nitrogen, from Fig. 1; X, nitrogen, from Fig. 4. 
B—oxygen, from Fig. 2. 


re Current inAmperes 


Fic. 5. Arc reignition potential for 
graphite in oxygen as a function of 
the arc current. The parameter is gas 
pressure in centimeters of mercury. 


ARCS 


Reignition Potential (volts) 


8 9 W 
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Fic. 6. Arc reignition potential for 
graphite in carbon dioxide as a 
function of arc current. The param- 
eter is gas pressure in centimeters 
of mercury. 


current at which the linear relationship between 
V, and pressure breaks down is indicated by the 
point at which B becomes zero. 

The curves of Fig. 8 show the effect of the 
initial rate of rise of impressed, or transformer 
voltage, upon the reignition potential. It is 
evident that the effect of increasing dV/dt in 


‘lowering V, is much less at low pressures than at 


the higher ones. The dotted portion of the curves 
is to indicate the sparking potential for the 
different pressures, i.e., the voltage approached 
as the rate of rise of impressed voltage following 
arc extinction is very small. The exact nature of 
the curve for the lower dV /di could not be de- 
termined. Although only a few points were 
possible, it is believed that the trend over this 
range is indicated. 


DISCUSSION 


The experimental results establish the fact that 
for the material used, and a circuit of small re- 
actance, the arc reignition potential is given, fora 
wide range of pressures, by V,=A+Bp. Since 
the field is not uniform, the curves of V, vs. 
pressure would be expected to depart somewhat 
from linearity at the higher pressures; the 
marked departure of carbon dioxide may be due 
to chemical effects. The results summarized in 
Table I apply, of course, only to the material and 
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Reignition Potential (volts) 


Initial Rate of Change of Voltage (Av /Sec.) 


Fic. 8. Are reignition potential of graphite in nitrogen 
as a function of the initial rate of rise of impressed voltage. 
The parameters are arc current and gas pressure. 


shape of electrodes used in this work; however 
they are of interest in that they indicate the 
variables present. The constant A follows the 
minimum sparking potential and quite probably 
the normal cathode drop, although this is not 
known for the materials and gas condition of the 
surfaces of the electrodes used. The fact that A is 
larger than the minimum sparking potential for 
oxygen and carbon dioxide may be due to gas 
impurities which are known to greatly affect the 
minimum sparking potential. A depends upon 
the rate of rise of impressed voltage, seeming to 
increase with this variable. This constant is 
independent of arc current except for carbon 
dioxide for which it increases as the current is 
decreased. The constant B depends on the nature 
of the gas, the current during the previous half- 
cycle of arcing, the frequency or dV dt, and quite 
probably on the gap and electrode material. It 
will be noted that the curves of B vs. J have the 
same form as the curves of I’, vs. J for the same 
gas. It is obvious that with the same gas and the 
same electrode materials this should be the case. 
Since the experimental results were quite uni- 
form, a knowledge of the value of the constant A 
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and of the equation for B makes it possible to 
predict the value of the arc reignition potential 
over a considerable range of current and pressure. 
However, since the constants of the B vs. I curve 
are purely empirical they give no clue to the par- 
ticular factors involved. 

The linear relation between V, and pressure, 
for the refractory electrodes used in nitrogen and 
in air, breaks down at a critical current for which 
the reignition potential is independent of pressure 
and has a value about equal to or somewhat 
greater than the normal cathode drop of a glow 
discharge. For currents larger than this value the 
reignition potential for graphite in air! decreases 
with increase in pressure; and the curves of Fig. 4 
indicate that this is probably true for nitrogen. 
The critical current indicates, probably, the 
approximate point at which the arc is reestab- 
lished without a preliminary glow; the residual 
thermionic emission due to heating by the previ- 
ous half-cycle of arcing being sufficient to restrike 
the arc with a voltage only a little above the arc 
burning voltage. 

For arc currents less than the critical value, a 
different situation exists. In this case thermionic 
emission, although of considerable magnitude for 
currents near the critical, operates in conjunction 
with a number of other factors in effecting re- 
ignition; all of which vary in relative importance 
with the state of the electrode surface, such as 
irregularities, impurities, adsorbed gases, etc. and 
also depend upon the nature of the gas, the 
pressure, etc. For the lower currents the reigni- 
tion potential for carbon, if pure, is approx- 
imately that of the relatively low melting point 
metals. Small amounts of low work function 
impurities will cause the reignition potential to 
vary from cycle to cycle, and larger amounts will 
lower it markedly.' The coefficient of heat con- 
ductivity, which probably varies with carbons of 
different manufacture, will also affect reignition 
by changing the rate of cooling of the surfaces 
after current zero. 

The effect of increasing the rate of rise of im- 
pressed voltage (increasing the transformer volt- 
age at constant frequency—assuming a sine 
wave) in lowering the reignition potential is 
probably due to the fact that less time is available 
for the cooling of the electrode surface. Also, less 
time is available for the deionizing of the gap. 
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Since the reactance of the transformer was 
fixed and the current adjusted by varying a series 
resistance, the power factor of the circuit ex- 
ternal to the arc itself, decreased as the r.m.s. 
value of arc current was increased. For the case 
where the effect of this constant circuit reactance 
would be greatest, i.e., for a large are current 
(variable resistance small) and low pressure 
(where the required voltage for reignition is 
small), it was found that the angle of lag was less 
than 6° and in this case reignition occurred 3.5° 
after zero current. Under these conditions ap- 
proximately 60 percent of the reignition potential 
was reached with a very high rate of rise of volt- 
age upon zero current and the balance at the 
relatively low rate of rise of the impressed volt- 
age.” For lower currents this effect was negligible 
because of the relatively high value of series 
resistance required. 

The results would seem to indicate that 
Paschen’s law can be extended to the reignition 
potential of the low pressure, short arc between 
refractory electrodes. It is necessary in this case, 
however, to specify the arc current since the 
constant B is different for each value of current. 
Recent glow discharge studies*: 4 have indicated 
that Paschen’s law holds for the reignition and 
burning potentials of the glow, and that the glow 
reignition potential is reduced when the rate of 
rise of applied voltage is increased.® Evidently the 
processes involved in the arc and glow reignitions 
are similar. 

It is probable that the early stages of the arc 
reignition are similar to the glow. discharge re- 
ignition. The difference between the two at the 
instant of zero current is due to the difference in 
the amount of ionization in the two cases. In the 
case of the arc a cathode spot must be formed. 
Hence the process of reignition must include a 
gaseous state in which conditions are such that 
this cathode spot may be formed. In previous 
papers': "it has been suggested that the reignition 
of an arc consists in the establishing of a ‘‘normal”’ 
glow discharge which quickly becomes ‘‘ab- 
normal,’’ the gap voltage and electrode current 
density increasing as the applied voltage in- 


* J. Slepian, Conduction of Electricity in Gases, Westing- 
house Course Number 38. 

*H. Fricke, Physik. Zeits. 34, 168 (1933). 

‘J. Thomson, Phil. Mag. 18, 696 (1934). 

* W. Spielhagen, Physik. Zeits. 34, 164 (1933). 

*S. S. Mackeown, Elect. Eng. 51, 386 (1932). 


creases. Then at some critical voltage the transi- 
tion to the arc is effected. This transition is 
known to take place very quickly if the condi- 
tions are favorable. The theory': ® proposes that 
the arc cathode spot forms as the result of the 
cumulative effects of positive ion bombardment, 
high electric field and high temperature on 
localized emission at the new cathode. Essen- 
tially the same suggestion has been made by 
Newman.’ 

The assumption of a preliminary abnormal 
glow provides a gas condition such that the 
transition voltage would be expected to vary just 
as it has been found to do in these experiments. 
First, there will be a considerable amount of 
positive ion bombardment and second, the posi- 
tive ion space charge of the glow will produce a 
high electric field at the cathode. The electron 
emission due to positive ion bombardment will 
depend upon the gas condition of the surface and 
the ion energy. The electrode surface being far 
from a true plane because of submicroscopic 
points, projecting crystal corners, etc., will be 
subjected to very high local electric fields. In 
these experiments the applied voltage will reach 
a value equal to the normal cathode drop in a 
time of one to three hundred microseconds (50 
cycles, 6600 and 2200 volts r.m.s.). When this 
point is reached the glow will start if the factors 
affecting the deionization of the gap have not 
been too active. These factors are recombination, 
diffusion, and reduction of thermionic emission of 
the electrodes due to cooling. 

It has been shown that in a few microseconds 
a short, turbulent, heavy current arc will recover 
a dielectric strength of several hundred volts per 
centimeter. The remainder of its dielectric 
strength is recovered at a very much lower rate. 
The results of Browne’s work are not directly 
applicable to the present case because they are 
for a longer arc, and an arc for which the turbu- 
lence and current were very much greater. How- 
ever, the effects of turbulence and current upon 
the dielectric strength are opposite to each other. 
Reducing the turbulence will reduce the di- 
electric strength of the gas at a given instant 
after zero current and reducing the current will 
increase the dielectric strength. Thus a rather 


*F.H. Newman, Phil. Mag. 18, 365 (1934). 
’ T. E. Browne, Jr., Physics 5, 103 (1934). 
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rough comparison of results may be made. It is 
evident that the dielectric strength indicated by 
the Browne results is exceeded early in the reigni- 
tion period in the present experiments, so that a 
glow discharge will result from the electric failure 
of the ionized gas. Further, it is interesting to 
note that the order of decreasing dielectric 
strength after one hundred microseconds is 
carbon dioxide, oxygen, and air, and that the 
order of decreasing reignition potential is the 
same. This is as it should be, since the gases hav- 
ing the highest recovered dielectric strength will 
require a higher voltage to initiate the pre- 
liminary discharge and will, no doubt, require a 
higher transition voltage as has been observed. 
The order indicated does not hold after thermi- 
onic emission has become a considerable factor 
since the experimental conditions have then been 
entirely altered. When the arc current is very 
small, it is necessary for the voltage to rise to the 
sparking value before the arc is struck. When the 
current is large, the thermionic emission for the 
graphite electrodes is high and the arc forms at a 
relatively low voltage without the preliminary 
glow. 

Recently published results would seem to 
justify the belief that a high electric field is an 
important factor in the glow-arc transition. The 
emission from a heated cathode may be in- 
creased by a factor of 10 by the action of a posi- 
tive ion space-charge field.’ It has been shown by 
Beams" that cold emission will start from impure 
mercury surfaces with average fields of 3.5 10° 
volts per centimeter and for pure mercury at 
1.8X 10° volts per centimeter; the latter value 
being considerably less than that required by the 
Fowler-Nordheim wave mechanics formula. It 
has been recently pointed out that emission at 
relatively low calculated fields could exist be- 
cause of points on the surface of the elec- 
trodes."': "8 In addition to the natural irregu- 
larities of the electrode surfaces, there is the 
effect produced by particles deposited upon the 
surface from the other electrode, or from other 
surfaces."* This effect will introduce an important 


*C.G. Found, Phys. Rev. 45, 526 (1934). 
*]. W. Beams, Phys. Rev. 44, 803 (1933). 
4. ]. Dempster, Phys. Rev. 46, 728 (1934). 
~, ©. Chambers, J. Frank. Inst. 218, 463 (1934). 
Tonks, Elect. Eng. 53, 239 (1934). 
H. Bennett, Phys. Rev. 44, 859 (1933). 
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change in the state of purity of the electrode 
surface as well as affecting the field, if the arc is 
studied at pressures below about 7 cm of mercury. 
Under this condition the arc cannot well be con- 
fined to the electrodes under study; but spreads 
to all exposed surfaces at different potentials and 
will therefore cause relatively large amounts of 
foreign material to be deposited upon the elec- 
trodes. Spectrograms recently taken at Harvard 
of the pure carbon before and after being arced at 
low pressures (5 cm of mercury) indicated that a 
considerable amount of copper and other im- 
purities had been deposited on the electrode 
surfaces, the nearest source being about twenty 
times the gap length from the gap. 

Another effect which cannot be overlooked is 
that which may result due to chemical reactions 
which may take place between the hot electrode 
vapor at and near the surface and the gas during 
the period of reignition. There is also the effect 
of dissociation of the gas and the effect of the 
ionization potential of any electrode vapor which 
may remain in the arc space after zero current. 
Just as all of these factors acting together affect 
the ionization or deionization of the gap, so will 
they lower or raise the voltage required to reignite 
the arc. A satisfactory theory of arc reignition 
must take into account all the factors outlined 
since they contribute in varying degrees to the 
reignition. 

The experimental results show the effect of 
increased current which increases the residual 
thermionic emission and the volume of ionized 
gas, and therefore lowers the reignition potential 
as required by the theory. A similar effect is due 
to increasing the rate of rise of impressed voltage. 
The close relation between the reignition poten- 
tial and the striking and burning voltages of the 
glow discharge also points to the correctness of 
the assumption of a preliminary abnormal glow. 
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When a dynamical system contains both mechanical and 
electromagnetic energies, the dynamical equations assume 
special forms. It is shown that when the dynamical system 
is nonholonomic, or when it is holonomic but is referred to 
nonholonomic reference axes, the equations still have the 
same form as those of a holonomic dynamical system and 
can be solved as such without any auxiliary equations. 
Following a suggestion of Lorentz, such dynamical systems 
are called here ‘‘quasi-holonomic’’ dynamical systems. 
Because of the special form of the equations, their develop- 
ment in tensor notation allows the introduction of concepts 
from the absolute calculus that apparently have not been 
hitherto applied in classical dynamics. It is shown here that 
the equation of motion of quasi-holonomic systems in 
tensor symbolism may be interpreted as specifying the 
trajectory of a point moving in an (n-dimensional) non- 


INTRODUCTION 


HE close analogy between the equations of 

performance of rotating electrical machin- 
ery and the equations of the differential geometry 
of abstract spaces (developed by the require- 
ments of modern electrodynamics) has been 
pointed out elsewhere,! illustrated with numerous 
examples. The key to the analogy was the 
recognition that the various interconnections of 
electrical networks and machines represent a 
group of transformations that leave several 
quadratic differential forms (stored magnetic 
energy, etc.) invariant. 

The performance of electromechanical systems 
may be represented by the motion of a point in 
an n-dimensional representative space. The fol- 
lowing types of spaces may be distinguished: 

(1) When the magnetic paths, the reference 
axes and the electric conductors are all sta- 
tionary, the space is a Euclidean space, having 
neither curvature nor torsion. 

(2) When the magnetic paths alone move, the 
space is a Riemannian space, having curvature, 
but no torsion. 

(3) When the reference axes alone move, the 
space is a space of distant parallelism, having 
torsion, but no curvature. 

(4) When the electric conductors alone move, 


_' Kron, J. Math. and Phys., 103-194 (1934); Kron, Gen. 
Elec. Rev. A serial beginning in the April, 1935 issue. 


Riemannian space, referred to holonomic or nonholonomic 
reference axes. When the field equations of Maxwell in 
tensor form are formulated for the electromagnetic part of 
the quasi-holonomic system, it is found that the electro- 
magnetic field must be expressed in terms of its components 
along an orthogonal ennuple, the field lying in a local 
Euclidean space, tangent to the underlying non-Riemannian 
space. By the use of such a nonholonomic reference frame the 
dynamical equations of Lagrange and the field equations of 
Maxwell are correlated. The equations show formal analogies 
with those of the five-dimensional unified field theory of 
relativity dynamics, where combined electromagnetic and 
gravitational systems are analyzed. Practical examples of 
quasi-holonomic systems are the numerous types of ro- 
tating electrical machinery used in industry. 


the space is a group space, having both curvature 
and torsion, but the components of the curva- 
ture and torsion tensors are all constants, satis- 
fying the fundamental structural relations of the 
theory of finite, continuous groups. 

(5) When the magnetic paths, the reference 
axes and the conductors all move with different 
velocities, the space is a general non-Riemannian 
space (metric space), with its curvature and 
torsion tensors satisfying the generalized identi- 
ties of Bianchi. 

When a cross section is made of any rotating 
machine during acceleration, at any one instant 
the current-density and flux-density waves are 
sinusoidally distributed in space. These space 
vectors may be considered as projections of 
vectors upon various two-dimensional planes, 
the vectors lying in a fangent n-dimensional 
Euclidean space and originating at the center of 
the rotor. That is, an electrical machine offers a 
visualization of vectors of the electromagnetic 
field lying in a centered local space, tangent to 
the underlying non-Riemannian (metric) space 
and moving along a trajectory in the underlying 
space. The translatory motion of the tangent 
space is accompanied by a rotation and slippage. 

In this paper the recently developed theory of 
nonholonomic reference frames? is used to put 


* Schouten and Struik, Einfiihrung in die neueren Metho- 
den der Differentialgeometrie (1935). The notation and ter- 
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the previously given dynamical equations in a 
more general form, so as to correlate them with 
the field equations of Maxwell, when they are 
applied to the electromagnetic part of a combined 
mechanical and electrical system. Just as in the 
unified field theory, no room can be found within 
the framework of a Riemannian space to represent 
the electromagnetic field in the presence of a 
mechanical sysiem. Hence even for classical 
dynamical systems the representative space must 
be non-Riemannian, to take care of the inter- 
actions between mechanical and electrical sys- 
tems. 


I. NONHOLONOMIC REFERENCE SYSTEMS 


(a) Let the geometric form of the equation of 
motion of a nonconservative holonomic dynamical 
system, expressed along holonomic reference 
axes, be given as 

d7x” 


-—-+[mn, k ] 
dt d 


dx” dx" 
——=f;,, (1) 
t dt 


Aim 


where a,» and f; are functions of x”. The holo- 
nomic form of the Christoffel symbol is defined as 


h 1 da ni OA mi OA mn 
Cn, #]~-( vias arenes ‘). (2) 


2\ ax" Ax" ax* 


Since frictional forces do not play any funda- 
mental part in the analysis, they will be omitted 
for the time being to make the equations simpler. 

(b) Let dx” be replaced by a new set of differ- 
entials dx* 


dx"=A,"dx* (3) 


by a nonholonomic relation, where A,” is a 
function of the old variables x” and where 


Ag" ox" ax, (4) 


That is the relation (3) is not integrable. As a 
consequence 


OA,,* Ax"HFOA,% OX”. (3) 


(c) By replacing also dim by @,aA,7A»” ete., 
the equation of motion expressed along the new 
nonholonomic reference system becomes 


d?x¢ | ; x dx8 f 
Aya— + LaB, y}— —=Jy 6) 
" di? Tr 


minology of this book is followed, except for the order of 
the indices of the various geometric objects which is kept 
the same as in Schouten’s Der Ricci-Kalkiil. 
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where the nonholonomic form of the Christoffel 
symbol is defined as 


[aB, y]=4(0a,3/0x°+0a,q/ dx8 — dq3/dX7) 
+246, at Qya, p—Qas, x 

where 2,3° is the ‘‘nonholonomic object”’ 
Q.3°=3Aq"Ag"(OA,°/ Ox" —OA ,*/ dx"), 


skew-symmetric in @ and 8, with a transforma- 
tion formula 


Qas? = 2,.°A QA ZA.! 


+A,°A3'(0A,°/axt—dAS dx), (9) 


and where 0a,3/0x* is written for (da,3/0x")Aq". 
In analogy to the Christoffel symbol the upper 
index is lowered as follows, 


(10) 


Qa, y= Qa3°057 


to simplify the form of the equations. 

In the n dynamical Eqs. (6) there are 2n 
unknowns; the n old variables x* and the 1 new 
differentials dx’, since the nonholonomic trans- 
formation (3) replaces only the old differentials, 
but not the old variables. The other » equations 
are the m equations of transformations (3). 


II. QuAsI-HOLONOMIC TRANSFORMATIONS 


(a) Let the m old variables x* be divided into 
two groups containing p and n—p variables and 
let the following assumption be made: 


(1) The components of the metric tensor a;,,, are functions 
only of the variables of the second group. 


Also let a group of transformations be intro- 
duced with the following additional assumptions: 


(2) The variables of the second group remain unchanged 
under the transformations to be considered, hence their 
differentials also remain unchanged. 

(3) The differentials of the variables of the first group are 
replaced by new differentials, in such a way that the 
elements in their matrix of transformation are functions 
only of the variables of the second group. 


(A more general case is to allow the variables of 
the second group to transform among themselves, 
instead of keeping them unchanged.) 

(b) Such a group of transformations occur in 
electrical machinery in going from reference axes 
moving with the conductors to axes moving with 
a different speed than the conductors. The first 
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ta 





Fic. 1. Rotating machines and their reference axes. 
a, holonomic, Riemannian ; 6, holonomic, non-Riemannian ; 
c, nonholonomic, non-Riemannian. Open circle = reference 
axis, connected to moving conductors; solid rectangle 
=reference axis, stationary in space; arrow =reference 
axis, moving at different speed than conductors. 


group of variables are the electric charges passing 
through the various conductors, the second group 
of variables are the displacements of : 


(1) The moving reference axes, 
(2) the moving rotor conductors, 
(3) the moving magnetic paths. 


A representative example is the transformation 
(Fig. 1, Section V) 
dx4+=dx?:, 


dx*=cos x*dx**+sin x*dx”, 
(11) 
dx’= —sin x*dx""+cos x*dx”, 
dx%=dx"%, 


dx*=dx’, 


with a matrix of transformation 
\a@ d, 


my 


d, o | 








a cos x* sin x* 





=h —sin x* cos x* 








Ss 0 


| 
| 

q. | 0 | 0 
| 


0 





where x* is the geometrical variable representing 
the rotor displacement. 

In general, quasi-holonomic transformations 
may occur in such dynamical systems in which 
both electrical and geometrical variables exist, 
since then dm,» is a function only of the geo- 
metrical variables. . 

Quasi-holonomic transformations are also used 
in the five-dimensional projective relativity’ 
where also: 


* Banesh Hoffman, Phys. Rev. 43, 615 (1933). 


DYNAMICAL 


SYSTEMS 145 


(1) The metric tensor is independent of the electrical 
coordinates, 

(2) in changing to a new reference system the space-time 
coordinates (geometrical variables) remain unchanged, 
or are transformed only among themselves, 

(3) a nonholonomic relation exists between the electrical 
axes, 

(4) the transformation matrix is independent of the 
electrical coordinates. 


(In the original five-dimensional unified field 
theory of Kaluza-Klein, the above cylindrical 
properties of the electrical variables also are 
assumed. However, the relation between the old 
and the new electrical variables is holonomic. 
For that reason the trajectory of the point in 
their theory is Riemannian. ) 

(c) After the transformation is made, the n 
new nonholonomic equations contain only n 
unknowns, the »—p variables of the second 
group and the p new differentials of the first 
group. In the new set of equations the old 
variables appear as new variables and the new 
equations can be solved as if they were holo- 
nomic. Hence, the new equations may be 
considered to describe the motion of a holonomic 
system. 


III. THE EQuaTION oF MOTION 


(a) Since the nonholonomic object Qas, , (and 
every other geometric object) is a function now 
only of the variables of the mew reference 
system, let its components be considered to form 
the components of another geometric object Saa,, 
a fensor of rank three, skew-symmetric in a and 8 


* 
Qas, Y = 


— Sapy- (13) 


In other words, let a new tensor S,3, be intro- 
duced whose components are identical (but 
negative) with the components of Q.s, , for this 
one reference system only, but different for other 
reference systems, since the transformation for- 
mula of Sag, is 


Sasy=Seewha'Ag’A ers (14) 


while that of Qa, , is the one given in Eq. (9). 

(b) Similarly, in this new reference system 
let [aB, y] of Eq. (7) be replaced by another 
geometric object I'as, , containing S,3, instead 
of Qas, y 


Pas, y= =[ap, ¥]+Sasy— Sypa— Syas, (15) 
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where now [a8, 7] is the holonomic Christoffel 
symbol defined in Eq. (2), since (da,3/0x")Aq” 
can be written in quasi-holonomic systems as 
0a ,8/ Ox". 

(c) If a new tensor is introduced to replace 
the three tensors 


T 087 = Sapy — Syba—Syas, (16) 
then the new affine connection is defined as 
Pas, = [ a, yj]+ T apy 


and the equation of motion (6) becomes for the 
new quasi-holonomic system :* 
d?x* dx* dx8 
Oya FT 09, oe OS 
dt? dt dt 


dyak*+[aB, y iti? + Task = f,. 


(17) 


(18) 


(19) 


(d) If the asymmetrical part of T'.3, , is taken 
in Eq. (15) then 


Pas}, y= Sasy, (20) 


since Sas, is skew-symmetric in its first two 
indices. Hence, the tensor S,.3, is identified as 
the torsion tensor used by Schouten and Struik. 
It also follows that the restriction 


Vals,y=0 (21) 


exists, since the affine connection in Eq. (15) 
does not contain the additional tensor Q,'%7 
=V,a*°". Hence the non-Riemannian space of 
quasi-holonomic systems is a me(ric space. 

(e) The effect of the frictional forces is to add 
an additional term R,.#*, so that the complete 
equation of motion of the quasi-holonomic 
system is 


Qyat*+Tag, xt +Ryat*=fy. (22) 


Such generalized trajectories have already 
been used in differential geometry.°® 

(f) In terms of fensors only, the equation of 
motion is 


647/5t+Ra1¢* = fr. (23) 


(g) In going over from this reference system 
to a third system by a quasi-holonomic transfor- 


‘ This equation with identical connection has been used 
for instance by Infeld, Zeits. f. Physik 50, 137 (1928), or 
by Banesh Hoffmann, reference 3, Eq. (13), etc. 

5 Veblen and Whitehead, ‘‘The Foundations of Differ- 
ential Geometry,’’ Cambridge Tracts in Math. and Math. 
Phys., No. 29, p. 72 (1932), Eq. 14.1. 
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mation, the equation of the new system is the 
same as that of the old, Eqs. (18) or (19), except 
now [a, y ] will be the nonholonomic Christoffel 
symbol defined in Eq. (7). Hence the affine 
connection I'as, , of this new system may be 
considered as being nonholonomic and non- 
Riemannian. 


IV. THE EQuaTIONS OF MOTION 
OF LAGRANGE 


(a) If, in the holonomic form of the La- 
grangian equation of motion 


d oT dT 


aS Fh 
dt 0%* odAx* 


T is replaced by (1/2)an,¢"%", the resulting 
equation is 
d?x™ (= 1 Od m n 


Qin—— + ss? as *) ane = fi. (24) 
di? Ox" 2 ax* : 


It should be noted that the term in the parenthe- 
sis is not the Christoffel symbol. However, it 
can be replaced by it since 


0a ml: 1 0a mk da nk 
——sn¢" =—-( —+ ans n, 
Ox" 2\ 0x" dx” 

(b) If a nonholonomic transformation is intro- 
duced in Eq. (24) instead of in Eq. (1), just as 
it was done in Section I, and a,» is replaced by 
AyaA,.IA»* etc., the resulting equation of motion 
has a different form from Eq. (6) 


| (== 1 ee) 
ax® 2 Ax? 
‘OA Pg OA ne 


+A rAar( naan 


Ox” Ox* 


ara ina =f, (25) 


If a,.%* is replaced by 07 /0%7, the resulting 
equation of motion is 


d ,oT oT 

= oS ae 

0A,° 
ox* 


oT 0A,2 
+—A,'A," ( 


0x? 


)e=F. (26) 
Ox" . 


This is the Boltzmann-Hamel form of the 
equation of motion valid for nonholonomic 
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reference axes as given by Whittaker, Schaeffer, 
etc.° 

(c) Eq. (25) can also be written for quasi- 
holonomic systems as 


d,at*+[oB, y]e2a®—2S,a3%%8=f,. (27) 


If a new tensor is introduced 


Pasy= —2Sya8 (28) 


skew-symmetric in its first and third indices, the 
equation of motion becomes 


Ayat*+LaB, y ]a*t®+ Pagyttt2=f,. (29) 


If a new connection is defined 


h 


h 
Aas, y= Lap, v¥]+Pasy _ Lap, 7] ei 2S ya3 (30) 
the dynamical form of the equaiion of motion is 


Ayat*+Aas, »t°¢? =f, (31) 


Both forms (31) and (18) of the equation of 
motion are given (without f,) in the unified field 
theory by the various writers as representing 
the trajectory of a charged particle. However 
they are not derived there from the equation of 
motion of Lagrange, but are set up from various 
geometrical, dynamical and philosophical con- 
siderations by a cut and try process. 

(d) Although the two forms of the equation of 
motion (18) and (31) have different affine 
connections (the geometric connection (15) and 
the dynamical connection (30)) both equations 
give identical answers, since from Eqs. (16) and 
(19) 

1. Sagyt%i8 =0, Sagy being skew-symmetric in its first two 
indices.’ 


2. (—Syag—Syga)t*iF = —2Syqgi%i* which is the same as 
Eq. (27). 
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However, in equations where the connections 
are not multiplied by z*z*, each of them appears 
to give a different answer. Such a case occurs 
for instance in trying to set up the field equations 
of Maxwell. 


(e) The tensor P,a3 is defined in Eq. (25) as 


Pass =A yAa"(dA 2 Ox"—0A a Ox" asa. ( 32) 


Let A,® be considered as a set of vectors A,. 
Since the rotation-coefficients of Ricci are defined 
as 


* R 
V7a8=A (y)*A (a)"V nA ass (33) 


(where the covariant derivative is taken with 
respect to the Christoffel symbol) the tensor 
Pya3 can be expressed in terms of them* as 


* 


P yap = ya8 — Yas: (34) 


Also if the Christoffel symbol is zero 


* 


Pas, y=Yasy (35) 


V. ROTATING ELECTRICAL MACHINERY 


(a) Among the numerous types of electrical 
machines the performances of those machines 
only are given by Eq. (1) (holonomic, Rieman- 
nian) where the reference axes are rigidly con- 
nected to the stationary or moving conductors. 
(Fig. 1a.) 

Its metric tensor da», represents all the self 
and mutual inductances of the various windings 
as a function of the rotor displacement angle x* 
and the moment of inertia of the rotor L. The 
resistance tensor R,,, (Eq. (22)) represents the 
resistances of the various windings. The metric 
tensor is 


b Is 





d, Mz cos x* 


— Mgsin x* 


0 





a M, cos x* 


La, cos? x*+L,, sin® x* 


(Lo,—La,) sin x* cos x* M, sin x* 





Gnn=b — My sin x* (Ly,—La,) sin x* cos x* 


M, cos x* 





qs 0 M, sin x" 





La, 





0 


Ss 0 














| 
| La, sin? x*+L,, cos? x* 
| 
| 
| 


M, cos x* | 
| 


0 0 








6 Whittaker, Analytical Dynamics, p. 43. ' 
* The fact that this term drops out of the equation of 
motion is also mentioned by Mandel, Zeits. f. Physik 56, 


838 (1929). 
5 McVittie, Phil. Mag. 8, 1033 (1930). 
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where d, and gq, are the stationary stator axes 
and a and b are the rotor reference axes moving 
with the same speed as the rotor. Axis s is the 
shaft axis along which all mechanical quantities 
(velocity, etc.) are represented. 

(b) Experience has shown that the equations 
of all industrial machinery come out in the 
simplest form if first the geometric objects of 
another machine are set up, whose rotor axes 
are stationary in space (Fig. 1b). This machine 
is called the ‘generalized machine’”’ and actually 
it contains several sets of windings on both 
stator and rotor. For this machine the salient 
pole is assumed to be stationary, although in 
the most general case it is moving. 

The transformation matrix A,” changing the 
geometric objects of machine (a) to that of 
machine (b) is given in Eq. (12). The metric ten- 
sor of the generalized machine is by d»,,Aq”A 3" 


SS 4 ode OG 
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the components of the tensor P.3, are (with s as 
a fixed index) 





























All the other components of P.s, are zero. (It is 
interesting to note that the Rotation-coefficients 
of Ricci do represent a rotation by changing the 
rotor flux-linkage wave to a rotor flux-density 
wave, at right angles to it in space.) 

(c) When a particular industrial machine is 
given (Fig. 2) its matrix of transformation A,* 
is read off the diagram of connections and its 
geometric objects, ad, Ro, [ oe, w] and P,.., are 
calculated from those of the “generalized ma- 
chine”’ given. 

When only one set of brushes exist on a rotor 
winding, the machine is a nonholonomic dy- 
namical system. 

(d) Calculation shows that in case of rotating 
electrical machinery: 


(1) Assuming first a holonomic machine (Fig. la, in which 
all reference axes move with the same speed as the 
conductors, but differ in speed from each other), in 
going over to a new nonholonomic reference frame 
there exists only one set of reference frames for which the 
nonholonomic object Qag,y (contained in [aB, y] of Eq. 
(6)) can be replaced by a tensor Pag,y. In that set of 
frames all axes move with identical speed. (Fig. 1b in 
which all reference axes are stationary.) 

The matrix of transformation A,* existing between the 
frames of this preferred set contain only constant com- 
ponents (the reference axes may form any constant 
angles with each other). 

In going over now from this preferred frame of reference 
to a general nonholonomic frame of reference, the non- 
holonomic object (contained in [aB, y] of Eq. (19)) can- 
not again be replaced by a tensor, since certain terms in 
the nonholonomic form of [a8, y] cancel. 


(3) 


Hence among all the possible sets of reference 
frames of rotating electrical machinery there are 
two preferred sets: 


(1) The reference axes are rigidly connected to the conductors. 
(2) The reference axes are rigidly connected to each other. 

These and other details will be covered in another 
publication. 
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VI. ADVANTAGES OF THE EQUATIONS IN 
TENSORIAL FORM 


In analyzing the performance of any given 
industrial machine with the aid of canonical 
equations, the engineer may use either the 
dynamical Eq. (26), or Eq. (29) in terms of 
tensor concepts. That is, the engineer may use 
either the holonomic machine (Fig. la) or the 
generalized machine (Fig. 1b) as his preferred 
reference machine. 

Experience has amply shown that the tensor form 
of the equation of motion is so far superior to the 
dynamical form for actual engineering calculations, 
that the use of the tensorial form is an absolute 
necessity in engineering and not jus! a convenience. 
All successful rotating-machine methods of anal- 
ysis (equivalent circuits, two-reaction theories, 
cross-field and revolving-field theories, etc.) have 
been based upon the tensorial form (Eq. (29)), 
though this fact is not generally realized. 

(1) In using the dynamical Eq. (26) holonomic 
(moving) reference axes n, k must always be 
assumed as the starting reference axes. That 
means that for machines with nonholonomic 
stationary axes (which comprises by far the 
largest number of industrial machines) the geo- 
metric objects of a machine with moving axes 
must be assumed first, containing always func- 
tions of @ (for instance Eq. (36)) and then they 
must be transformed with the aid of a transfor- 
mation matrix A,” (for instance Eq. (12)) whose 
components are also functions of 6. 

In other words, the dynamical Eq. (26) compel 
the engineer to select as his starting reference 
system one, whose equations happen to assume 
a most complicated form. The change from this 
complex form to the final simple form represents 
extra mental and physical work for the engineer. 

On the other hand the tensorial form of the 
equations allows the engineer to select as his 
starting (preferred) reference system any one in 
which the equations assume the simplest form. 
That reference system happens to be a quasi- 
holonomic system, that of the “generalized 
machine” (Fig. 1b) whose geometric objects are 
of very simple form and are not functions of 86, 
as shown in the previous section, Eqs. (37), (38) 
and (39), 


The transformation matrix A,* representing 
the step from this generalized machine to an 
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Fic. 2. Repulsion motor and its matrix of transformation. 


actual industrial machine (Fig. 2) (that is to 
any allowable reference system) is of the simplest 
form, usually containing only constant compo- 
nents, hence the calculation of the geometric 
objects of the actual machine is of simple 
character. 

(2) In using the dynamical Eq. (26) the 
matrix of transformation A,” in many Cases is a 
singular matrix and the inverse transformation 
A»*, which is needed in (26), cannot be estab- 
lished. In such cases nonholonomic ‘‘equaiions of 
constraints’’ must be set up in addition to the 
holonomic dynamical equations and the calcu- 
lations become still more involved. 

However, in transforming the tensorial Eq. 
(29) of the generalized machine, the inverse of 
a singular matrix of transformation is not needed, 
since the indices of all geometric objects are 
covariant indices. 

(3) The physical picture of the phenomena 
taking place inside the rotating machine, as 
given by the various geometric objects and 
terms, depends on the starting reference system. 
When a holonomic machine is assumed as the 
starting reference system, the physical interpre- 
tation of each term for any other machine is 
most unusual, artificial and is of no aid to the 
engineer in his visualization. 

On the other hand the physical picture given by 
the various terms of the tensorial Eq. (29) is the 
same that has been used by electrical engineers 
previously in setting up ihe equations of a machine 
from purely physical considerations, without the 
use of any canonical equations. 

(4) It will be shown later in Section VII that 
the tensor form of Maxwell’s field equations also 
makes it necessary to consider the extra term in 
Eq. (26) as a tensor P.3,<*%%* instead of a non- 
tensor term 20,3, ,#°%*, in order to correlate the 
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dynamical equations of Lagrange with the field 
equations of Maxwell. 

There are still other considerations, all pointing 
to the advantage of the tensorial formulation of 
engineering problems. 


VIL. Tue Fretp Equations of MAXWELL 


(a) If a cross section is made of any rotating 
machine, it is found that the current-density and 
flux-density waves are sinusoidally distributed in 
space (neglecting higher space harmonics). This 
suggests trying to set up Maxwell's field equa- 
tions in tensor form for the electromagnetic part 
of the quasi-holonomic system under consider- 
ation. It is of course known that Maxwell's 
equations do not incorporate the effect of 
mechanical energies and are insufficient to deter- 
mine the performance of electrical machinery 
during acceleration. 

Calculation shows that the field equations of 
Maxwell for a quasi-holonomic system can be 
set up only when the tensorial equation of motion 
of Lagrange is used, Eqs. (19) or (29) which is 
derived by assuming a special nonholonomic 
reference frame as a starting point, in which all 
axes have the same velocity. 

It is found that in order to get correct answers, 
both connections I'as, , and Ags, » have to be 
considered. Covariant differentiation with re- 


G 
spect to I'as, , will be denoted by V (G=geomet- 
rical connection) and that with respect to Ags, , 
D 
will be denoted by V (D =dynamical connection). 


Since only the electromagnetic part of the system 
is considered, the moment of inertia of the rotor 
L is replaced by unity. 

It should be noted that since the electrostatic 
field existing in the dynamical system is not 
considered, the time axis does not appear in the 
analysis. 

(b) The flux-linkage wave (electrokinetic mo- 
mentum) of all axes is defined by Maxwell’ as the 
electromagnetic vector-potential. However it is 
necessary to assume, that this definition is valid 
only for the generalized machine with stationary 


axes. That is 
Ya = desi? (40) 


® Maxwell, Electricity and Magnetism, Vol. 2, p. 247. 


GABRIEL KRON 


For the non-Riemannian generalized machine 
(Fig. 1b). 
d, d, qr qs 


wt | j 
a= ad, 





Yad, | a, | 4, 


| # | (41) 





where gu,=14%Ma+i"La,, also ¢4,=i1%M,+i%L,, 
etc. 

For the Riemannian machine (Fig. 1a) ¢,, 
must be defined in terms of ¢,, instead of a,,,i” 


Om = Yad m% (42) 


where A,” is given in Eq. (12). 
d; a b 





¢a, COS X*+ —¢a, sin x* | 
=> | 
| ¢q, Sin x* + yy, cos x* | 


| 





(c) The first equation of Maxwell is 


G 
Fag =2V (as). 
Since the curl of a vector is 


G h 
V (a¥3} = 9[a¥3}] — Daa ?¢y, 


(45) 
the skew-symmetric electromagnetic tensor is! 


Fag=Oga/ 0x8 —Ayg/dx*—2Si37%¢,. (46) 


For the non-Riemannian generalized machine 
(Fig. 1b), the components of ¢, being constants, 


(47) 


Fi3= ae 28037 Py _ P ps0k*. 
Calculation shows that the components of the 
skew-symmetric field tensor are 
\B 

















representing (twice) the resultant flux-densily 
wave of the rotor, at right angles in space to the 
flux-linkage wave. 

Similar formula for Fag is derived from a variational 


principle by Novobatzky, Zeits. f. Physik 89, 379 (1934), 
Eq. (25). 
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It should be noted that if in Eq. (45) the 
tensor — 53,37 had been replaced by its original 
definition Qas7 (Eq. (13)), a nontensor object, 
then the right-hand side of Eq. (47) could not 
have been defined as F,3, a tensor. Hence, the re- 
placement of Qag? by — Sig’ is required by Max- 
well’s equation. 

For the Riemannian machine (Fig. 1a) S337 is 
zero and 


Fruin=O¢m/Ox" —OGn Ox”. (49) 
By differentiating Eq. (43) 
m» 


d. 








a 


Fun n = b 








where 


a= $a, COS X*+ vy, Sin x* 
(51) 
¢b= — ga, SIN X*+ Gy, COS X* 


giving the rotor flux-density wave. F,,, could 
have also been found from Fy; (Eq. (48)) by a 
tensor transformation F,,,= FasAm*A,®, where 
A,” is given in Eq. (12). 

(d) If Eq. (47) is substituted into the dy- 
namical equation of motion (29), it becomes for 
any machine or quasi-holonomic system 


Ayakt*+[aB, y ]a*e?+ Fyat* = fy. (52) 


This form of the equation of motion of a charged 
particle is given (without f,) in many papers on 
the unified field theory." It should be noted that 
in the Riemannian equation of motion (1) there 
is no place left for the electromagnetic-field 
quantities. 

In case of rotating machinery the term F,,%* 
represents: 


(1) The total torque exerted on the rotor (y=s) 
(2) The rotor generated voltages due to rotation of the 
conductors. 


"See the paper by Mandel, reference 7, Eqs. (7) and 
(10) or Tolman, Relativity, Thermodynamics and Cosmology 
(1935), Eq. (103.1), ete. 
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The term [a6,y ]é*z* gives the voltages gen- 
erated in all stator and rotor windings due to 
the motion of the reference axes and of the 
magnetic paths. It also gives the torque on the 
latter. 
(e) From Eq. (44) it automatically follows that 


G 
ViaFsy =0. (53) 


(f) The second independent equation of Maxwell 
is 
D 


V"F.s=Ss ( 4) 
53=0F*%B/Ax*+Ag,, FY. (55) 
For the non-Riemannian generalized machine 


S3=Agy, fv*=Ps,,F. (56) 


Calculation shows that 


ds d,; qr qs 


| Nl 
| 9 | ea, | @, 0 





5p = | 





representing the covariant rolor currenis in the 
stationary rotor axes. (See Eq. (40).) 
For the Riemannian machine the last term of 
Eq. (55) is zero and 
Su = OF*,,,/Ox". (58) 
Differentiating Eqs. (50) and (51) 
ds da b qs s 


«| 0 0 | (59) 





Sm = | 0 Ya 





representing the covariant rotor currenis in the 
moving rotor axes. 

(g) From Eq. (52) follows the force equation 
(torque on the rotor) 


Fgs* = Da. (60) 


(h) From Eq. (54) it automatically follows 
that 
D 


V%S,=0. (61) 


(i) It should be noted that Maxwell's equa- 
tions do not define geometric objects for the 
staior axes and F,g and s, have no components 
along d, and q,. That is, the free indices in 
Maxwell’s equations do not assume stator values 
(but the dummy indices do). 
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(j) The stress-energy tensor is defined as 


Ta3= F**F 3,—\ga8F oF. (62) 


Defining ga3= >. mAa"A 3” for the rotor axes only, 
(Eq. (12)) it gives the unit matrix 6a.” 

For the non-Riemannian machine calculation 
gives the symmetrical tensor 


d, d, 





0 | —~ Pd,Par 

















where E is the magnetic energy stored in the 


rotor 
; gar gaye 
E=- + 
ag | 


and Ly,’ and L,,’ are the inductances of the rotor 
windings measured with the stator windings 
short-circuited 


L,’ =(L.L,— M’*)/L,. 


(64) 


(65) 
The stress-energy tensor represents: 


2 Levi-Civita, A Simplified Presentation of Einstein's 
Unified Field Equations (Blackie & Son, 1929), Eq. (35). 


GABRIEL 


KRON 


(1) The Maxwell stresses existing between the rotor 
windings. 

(2) The magnetic energy stored in the rotor. 

(k) The existence of two (and only two) sets of 
preferred reference systems brings into sharp con- 
trast the two schools of thought that exist in 
analyzing rotating machine performance: namely 
those using flux linkages and those using flux 
densities: 

(1) In starting with reference axes moving with the con- 

ductors, the machine analysis becomes a problem in 
particle dynamics, using the holonomic equation of 
motion of Lagrange, Eq. (1) as a starting point. In the 
analysis only quantities measured at the machine 
terminals occur, the flux-density wave, the current- 
density wave and in general the field equations of 
Maxwell are nonexistent. 
In starting with stationary reference axes, moving 
sinusoidal waves inside the machine appear, that were 
nonexistent in the previous analysis. The quasi- 
holonomic equations of motion, Eqs. (52) or (19) are 
used as a Starting point and the field equations of 
Maxwell apply to the space waves. 


Hence the dual point of view, that is, considera- 
tion of a rotating electrical machine as a collec- 
tion of particles (terminal quantities) or as a 
collection of waves (sinusoidal space quantities), 
is a consequence of the existence of two preferred 
reference frames from which the phenomenon 
may be viewed. 

The holonomic preferred frame must be used as 
a starting point for the dynamical equations of 
Lagrange, the quasi-holonomic preferred frame 
must be used as a starting point for the field 
equations of Maxwell. 





PHYSICS 


VOLUME 7 


The Thermal Conductivity of Tungsten 


Water C. MICHELS AND Martua C ox, Department of Physics, Bryn Mawr College 
(Received December 28, 1935) 


The theory of the heat losses from an electrically heated wire, as given by Roberts and 
Kannuluik, has been extended to include a wire supported by springs. The results obtained have 
been used in a measurement of the thermal conductivity of commercial tungsten wire in the 
temperature range from 78°K to 273°K. There are definite indications that the Wiedemann- 


Franz law breaks down badly in this region. 





N the course of an investigation of accommoda- 

tion coefficients, recently carried out in this 
laboratory,! it became necessary to use the ther- 
mal conductivity of the tungsten wires in order to 
correct for the end losses in the Pirani gauge used. 
A search of the literature revealed that, although 
tungsten has been very thoroughly studied at 
high temperatures,’ very little work has been 
done on its thermal properties in the neighbor- 
- hood of room temperature. Kannuluik® recently 
reported the value of the thermal conductivity at 
273°K as 1.67 watts/cm deg., checking quite well 
with the older values, and Roberts,‘ in the course 
of his work on accommodation coefficients, 
estimated the value at 80°K to be 2.6 watts/cm 
deg. No determinations were found at other 
temperatures in this neighborhood. Roberts’ 
value, moreover, failed to check with a rough 
determination made in the course of our accom- 
modation coefficient work. It became necessary, 
therefore, to make a set of determinations of the 
conductivity, and we chose the measurement of 
the losses from an electrically heated wire in 
vacuum as the most convenient and reliable 
method. The values which are reported in this 
article are essentially preliminary, and it is 
planned to carry out further determinations over 
a wider temperature range and with improved 
apparatus, but we feel justified in presenting 
them because of the paucity of existing data. 


I. THEORY OF AN ELECTRICALLY HEATED 
WIRE WITH END SUPPORTS 


Kannuluik*® and Roberts‘ have developed the 
theory of an electrically heated wire with its 


' Michels and White, Phys. Rev. 47, 197A (1935). 
wat” see Forsythe and Worthing, Astrophys. J. 61, 146 
{ ). 
* Kannuluik, Proc. Roy. Soc. A131, 320 (1931). 
* Roberts, Proc. Roy. Soc. A135, 192 (1932). 


ends and surroundings at a common temperature, 
T. In practice, however, the heat treatment of 
the wire requires that it be supported by springs,® 
and the temperature drop in these introduces an 
appreciable error in some cases. 

Consider a uniform cylindrical wire of radius a, 
length / and resistivity p, mounted in vacuum, 
the surroundings being at a temperature 7. Sup- 
pose that the wire is supported by two springs 
made of the same material, each of length /, and 
radius a,. Finally, let the current flowing in the 
wire be z and the temperature drop in the springs 
(which is to be evaluated later) be AT. In a sec- 
tion of the wire of length dx, the electrical power 
input is 7pdx/ma?. If the temperature excess 
above the ends is ¢, the net radiation loss is 


2raeo | (T+1+AT)*—T* I} dx, 


where ¢ is the emissivity of the wire and ga is the 
Stefan-Boltzmann constant. When (t+ AT7)/T is 
small, this reduces to 


8raeoT*(t+ AT )dx. 
Finally, the conduction loss from the section is 
equal to 

— na’K (d*t/dx*)dx, 
where K is the thermal conductivity of the wire. 


Equating the power input to the sum of these 
two losses, we have: 


d*t ‘dx? -—At+B=0, (1) 
where A=8eo7*/aK and B=7*p/r’atK—AAT. 


Solving this equation, subject to the boundary 
conditions that ‘=0 at x=0 and at x=1, we 
obtain the temperature distribution of the wire: 


sinh xA!+sinh (J—x)A}) 


t=— [° (2) 


A sinh /A} 
® Blodgett and Langmuir, Rev. Sci. Inst. 5, 321 (1934). 
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This expression for ¢ enables us to find the total 
conduction and radiation loss from the wire. 
First, setting the conduction loss equal to 


2na*K (di /dx) ,~», 
O,.=2na*KBf(A)/A}, (3) 


where f(A) has the value (cosh /A}—1) /sinh /A}. 
Finally, the total radiation loss is: 


we find: 


QO, =8naeoT*(t+ AT)I, (4) 


where / is the average value of ¢, obtained by 
integrating Eq. (2) over the length of the wire. 

It is interesting to note that Eq. (2) gives a 
parabolic temperature distribution if the radia- 
tion or lateral losses are negligible, and a rec- 
tangular distribution if the conduction losses 
are negligible. For the former the temperature 
gradient at the ends has a magnitude B/ /2, while 
for the latter it is infinite. In the case of any 
actual wire the value of this gradient will be 
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between these extremes, and will increase as the 
lateral losses are increased. Consequently, if the 
lateral losses are made greater by any means, 
such as increasing the emissivity or surrounding 
the wire with gas, the conduction losses also 
increase for a given average temperature rise. This 
is an extremely important factor in determining 
the accuracy of any measurement involving the 
heat losses from a wire. 

We shall now consider the springs themselves, 
and shall use the subscript s to denote quantities 
referring to the springs, otherwise keeping our 
earlier notation. If we let 


A,=8e,0T*/a,K, and B,=i*p,/7r’a3K, 


Eq. (1) represents the springs. The boundary 
conditions are now /,=0 when x=0, and 
Q,./2=7a/7Kdi./dx at x=l,. 


Imposing these conditions, we find that: 


(B,/A,)(cosh/,4,!—1)—Q, 27a2KA.} sinh 1,4.) 


AT = (t.) 2=1,=—————_ 


(B.l,/A,) sinh1,A 3—(B,/ A.) coshl,A,)+ (QO, /27a,K;A,') (cosh 1,A .—1) 


‘= ——S 


In case either the emissivity or the temperature 
low that the radiation may be 
neglected, these equations become simplified. In 
this case f(A)/A! approaches the value / 2 and 
we have: 


is so losses 


(=(B/2)(lx—x°*), 

i= BF /2, 
AT=BJd2/2+Qil./20a7K, 

(,=Bd,4+Q.. 4racK. 


and 


II. DETERMINATION OF THE THERMAL 
CONDUCTIVITY 


The gauge used for these measurements con- 
sisted of a tungsten wire of 0.00254 cm radius and 
14.8 cm long mounted in a Pyrex tube, between 
springs of 0.00635 cm radius wire. Each spring 
contained 6.0 cm of wire. The tube was so con- 


/, cosh 1,4. 


cosh /,A a 





(6) 





structed that it could be lowered into a one liter 
Dewar flask for thermostatic purposes. After a 
thorough baking on the pumps (mercury pumps 
with CO, and acetone as a freezing agent), the 
wire was aged by flashing the filament at white 
heat for several hours. This flashing removes the 
tungsten oxide from the surface and etches the 
surface to a stable condition.':* The tube was 
then kept on the pumps at a measured pressure 
of <10-* mm during the measurements. 

A standard resistance of one ohm was con- 
nected in series with the tube and the potential 
differences across this resistance and across the 
tube itself were measured by a potentiometer. A 
series of currents were used, the maximum giving 
an increase of resistance of about 25 percent. At 
two of the three temperatures used, those ob- 
tained by thermostating with melting ice and a 


6 Taylor and Langmuir, Phys. Rev. 44, 423 (1933). 
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carbon dioxide snow-acetone mixture, the resist- 
ance of the tube was found to be a linear function 
of the power input, but the departure from 
linearity at the third temperature (that of liquid 
nitrogen) was distinctly outside the experimental 
error. Since the resistance is very closely a linear 
function of temperature and the power input is 
the sum of Q, and Q,, we should expect from 
Eqs. (3) and (4) that this linear relation would 
hold. Two approximations, however, have been 
made in deriving these equations, when it was 
assumed that the resistivity is a constant over the 
length of the wire and that the radiation loss is 
proportional to the temperature excess above the 
surroundings. These two errors compensate 
somewhat, since the increase in resistivity ac- 
counts for a higher power input to those sections 
of the filament from which the radiation loss is 
greater than we have assumed. But at the 
temperature of liquid nitrogen the radiation loss 
is negligible, and the departure from linearity is 
undoubtedly due to variation of resistivity with 
temperature. 

From the measured resistance and power input 
the thermal conductivity was found by successive 
approximations. Assuming a value of the thermal 
conductivity and taking the emissivity of the 
tungsten oxide on the springs as ten times that of 
clean tungsten, a consistent solution of Eq. (6) 
was found. From this the resistance of the 
springs could be determined, so that the resist- 
ance and the average temperature rise of the 
filament itself were obtainable. Hence Q, could 
be calculated from Eq. (4) and subtracted from 
the power input to the filament to give Q,. 
Putting this quantity into Eq. (5), AT was 
found. Putting this into Eq. (3), f(A)/A! was 
obtained and the value of K was found graphi- 
cally from a plot against this quantity. This 
method was repeated until the assumed value of 
K and the final value agreed. Values of the 
emissivity were taken from the data of Langmuir 
and Taylor’ and the resistivity from those of 


* Langmuir and Taylor, J. Opt. Soc. Am. 25, 321 (1935). 
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TABLE I. 











K/oT 
(e.s.u.?) 
2.36 K 107" 
2.98 
3.87 


TEMPERATURE 
(deg. K) 


THERMAL CONDUCTIVITY 
(watts/cm deg.) 


78 2.26 
194 1.60 
273 1.88 








Holburn,*® with which our own measurements 
agreed after a small correction was made for the 
resistances of the joints. At the lowest tempera- 
ture Eqs. (8)—(10) were used for the calculation. 

The values of the thermal conductivity found 
in this way are given in Table I. Because of the 
large influence of the diameter of the wire on the 
value of K, we estimate the uncertainty in these 
values to be about ten percent. In the third 
column of the table, the values of the Wiede- 
mann-Franz “‘constant’”’ are tabulated. In con- 
nection with this quantity it is interesting to note 
that the value which we obtain for it checks with 
that obtained by Kannaluik much more closely 
than do the values of the thermal conductivities 
obtained in the two experiments, indicating that 
its value is not seriously affected by impurities 
which, apparently, change the electrical and 
thermal conductivities in the same ratio, as 
would be expected from theoretical considera- 
tions. A very serious departure from the ordinary 
theory of the Wiedemann-Franz law is indicated, 
however, in that the ratio K/oT seems to be 
approaching a value of zero as the temperature 
approaches zero. 

Further work on the thermal conductivities of 
the metals in this temperature range is planned 
with two objects in view, to obtain more precise 
values for practical application and to test this 
interesting and peculiar behavior of the Wiede- 
mann-Franz constant. We expect to have avail- 
able for this work a sample of extremely pure 
tungsten wire. 

The authors wish to express their gratitude for 
a gift of the late Marion E. Reilly which made 
possible the purchase of some of the apparatus 
used in this work. 


8 Holburn, International Critical Tables, Vol. 6, p. 129. 
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The Thermal Stress in a Strip Due to Variation of Temperature Along the Length 
and Through the Thickness 


J. N. Goopier, Ontario Reseach Foundation, Toronto 
(Received February 6, 1936) 


The thermal stress developed by unequal heating of a 
flat strip is investigated, the temperature varying along the 
length and through the thickness, but not across the width. 
The problem is reduced to one of edge tractions only, for 
which general solutions are already known. From the 
properties of these solutions it is shown that for heating 
local to a narrow transverse band across the width, as might 
occur in welding processes, there is a concentration of stress 
near the edge of the strip, at the hottest part, and a 
maximum tensile stress of magnitude EaTjyax,, when the 


temperature is constant through the thickness. The 
corresponding quantity is found also for temperatures vary- 
ing through the thickness. For steel EaTnax, =350Tinax. 
pounds per in.?, the temperature being measured in centi- 
grade degrees above that of the cool parts of the strip, and 
this approaches the yield point of mild steel for a maximum 
temperature (rise) of the order of 150°C. For a given 
maximum temperature, the maximum stress is reduced by 
about 40 percent by keeping one face cold. 





1. 


HEN a solid is unevenly heated, the 
tendency of the hotter parts to expand is 
restrained by the rigidity of the cooler parts, and 

a state of thermal stress is developed. 

The particular case of the thin flat strip heated 
locally all across the width, is suggested by weld- 
ing processes. Two types of temperature distribu- 
tion will be considered. 

(i) The temperature, 7, is constant through the thickness, 
and across the width. With axes as in Fig. 1, it isa 
function of x only. 

(ii) The temperature on either face (= +}h) is a function 
of x only, and the variation through the thickness / 
from 7\(x) to 72(x) is linear. The temperature at any 
point can then be expressed as 


T= }(T.4+T2)+(2/h)(TM— T2). 


This type of distribution will occur when one face of the 
strip only is heated. 
° | 


b 


| 














Fic. 1. 
2. First Type OF TEMPERATURE DISTRIBUTION 


The distribution (i) will give rise to a state of 
generalized plane stress. Let u, v be the displace- 


ments, X,, Y,, X, stress components, averaged 
across the thickness. The stress-strain relations 
are 


E fou Ov ' 
—+o—-—(1+o)aT 
1—o’ Lox ov J 


X,= 


E foav_ ou 7 
Y,=——— | —+o——(1+o0)aT 
1—o*Ldvy Ox 


E ‘Ov Ou 
ng ie 
2(1+¢)\dx dy 


where EF is Young’s modulus, o Poisson’s ratio, 
and a@ the coefficient of linear expansion. When 
these are substituted in the equations of equi- 
librium, the result is the two equations 


ne) Ou dv l-co/? @& 
SANDE tan 
dx dy/ \dx dy 1+o\dx? dy 








When 7 is a function of x only, we find a particu- 
lar integral by assuming that wu is a function of x 
only, and v=0. This leads to du/dx=(1+0)aT, 
and to the stress components 


X,=X,=0, Y,=—EaT. (1) 


There is therefore normal traction —EaT on the 
edges y=0, b. In order to make these edges free 
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Fic. 2. 


it is necessary to superpose on (1) the stress distri- 
bution due to equal and opposite edge tractions 
(Y,=EaT on y=0, b) and the resultant distribu- 
tion will then be the thermal stress in the strip. 

The problem of the strip of finite width loaded 
by edge forces has been completely solved,! and 
formulas for the thermal stress can now be 
written down. 


3. 


We shall consider in detail only the case when 
the region of high temperature is narrow, in the x 
direction, in comparison with the width of the 
strip. The traction EaT is then concentrated on 
a small length of the edge, and the stress it 
produces will be local to this length. The exist- 
ence of the other edge is immaterial, and each 
half of the strip can be treated as a half-plane. 
The stress in the half-plane due to normal 
traction p constant over a segment AB of the 
edge y=0 (Fig. 2) is derivable® from the Airy 
stress function x = (p/27)(r:?@:—72°82) by means 
of the formulas 


K,=0x/dy’, Y,=0"x/dx*, X,=—0x/dxdy. 


Then 


X,=(p 27)[2(@;—62)+sin 26,;—sin 262 |, 
Y,=(p/27)[2(0:—62) —sin 26,+sin 262 |, (2) 
X,=(p/27)[cos 26;—cos 26. ]. 


It is clear that these satisfy the boundary condi- 
tion: Y,=p within AB, zero outside AB; and 
X,=0. 

Making the segment infinitesimal, and in- 
troducing ~ as the coordinate of A measured 


'L. N. G. Filon, Phil. Trans. Roy. Soc. A201, 63 (1903); 
and T. von Karman, Abh. Aerodynam. Inst., Tech. Hoch- 
schule, Aachen 7, 3 (1927). 

* J. H. Michell, Proc. Lond. Math. Soc. 34, 139 (1901). 
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along the edge, the stress function becomes 
(p/2m)-(0/d€)r6’-dé where r” =(x—£)?+ y?, and 
6’ is the principal value of tan~! (y/x—&). Carry- 
ing out the differentiation, and integrating the 
effects of all the elements of a load function p(£) 
we find that the corresponding stress function is 


’ 


1 +o 
x(x, y)= --f p(£)(x—£) tan-' ——dé (3) 
Te _.., 


x—& 


except for a part which makes no contribution to 
the stress. Traction EaT is obtained by sub- 
stituting p(t) = EaT(é). 

From (2) we obtain, for the uniformly loaded 
segment, X,=Y, on the edge, where Y, is of 
course equal to the applied traction p. It follows 
that at the edge X, is zero outside the loaded 
segment; when we assemble different loadings on 
a series of infinitesimal segments, so as to build 
up the arbitrary loading p(£) or p(x), the stress 
component X,, at the edge, within each infini- 
tesimal segment is simply due to the load on the 
segment itself, and we shall have in the resultant 
stress distribution, defined by (3), (Xz),-0= p(x). 

The normal edge traction EaT therefore pro- 
duces an edge stress X,=EaT. There is no 
contribution to X, from (1), so that this is the 
actual component of the thermal stress. The 
heating being local, the stress due to the edge 
tractions diminishes rapidly from the edge in- 
wards. If the temperature rise is appreciable only 
between x= +) the stress-component X, will be 
negligible, in comparison with the edge value, 
beyond distances from the edge of the order 2. 
In the interior the principal part of the stress will 
be given by (1). 

The character of the components X, and Y, 
across the middle section (x=0) is indicated in 
Fig. 3. The component X, is tensile at the edges 
with a maximum value EaTmax.. For steel this 
is about 3507 ax. pounds per in.?, Tinax. being 
measured in degrees centigrade above the un- 
heated part of the strip. The degree of heating 
required for a tensile stress of 50,000 Ibs. per sq. 
in.—near the yield point of an ordinary steel, is 
given by Tinax.=143°C. Further heating would 
lead, in such a material, to permanent deforma- 
tion, and a residual stress would remain after 
cooling. 
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4. SECOND TYPE OF TEMPERATURE DISTRIBUTION 


The second type of temperature distribution, 
(ii) of $1, involves both expansion of the strip in 
its own plane, and flexure transverse to it. The 
expansion, associated with the mean temperature 
3(7,:+T:2), produces stress of the type already 
discussed, and is independent of the flexure. But 
the flexure, associated with the temperature 
difference, 7,—7>, will in general be influenced 
by the mean stress due to the mean temperature. 

To form the differential equations of the flex- 
ure, let G,, G. be the bending stress couples per 
unit length on sections normal to the x and y 
axes, respectively, /7,, the torsional stress couple 
on a section normal to the x axis, N;, Ne the 
shearing stress resultants per unit length, and D 
the flexural rigidity Eh’ /12(1—o?). The equations 
of equilibrium are 


0G; ‘Ox — oll, oy _ N, = 0, 
0Gs ‘Oy —OH,/dx— N2=0, (3) 


ON 1 0 Ne rw rw 
+X,-—+2X,——_-+ Y,—-=0, ' 
Ox dy Ox? Oxdy P 

the last expressing the equilibrium of the trans- 
verse force components (i.e., in the z direction) ; 
X,, X,, and Y, are the mean stress components 
due to the mean temperature, and w is the 
transverse deflection. The relations between the 
stress couples, stress resultants, and the deflec- 


tion w are 


Pw fw a i 
G,= — D| —-+o0—- —-(1+¢)(7T,—Th?2) 
Ldx> ay’ oh 


ew Pw a 
—+o——-(1+0)(Ti—T>) |, 
Ldy> ax?” sh 


Ge= —D 








4 


rw 
IT, = D(1 —«)——, 
axdy 


. ‘ 
v2w—(1+0)(7,—T2) |, 


dx L h 


- 4 


a 
v?w——(1+0)(7,—T2) |. 
h 4 








The formulas for G,, Ge, for instance, are derived 
from the consideration that the unrestrained 
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Fic. 3. Character of stress distribution in hottest section 
of locally heated strip. 


thermal expansion of an element of the strip 
would produce spherical curvature (7,—72)a h. 
The couples are responsible not for the whole 
actual curvature, but only for the difference 
between it and the free thermal curvature. 
Elimination of Gi, Ge, H1, Ni, Ne, from (3) and 
(4) leads to 


rw 
D-viw—X,—--2 


Ox? Oxdy oy" 


a 
wi +o)¥7(1,;—T2). (5) 
i 


vw stands for e* ax?-+02/dy? and viev’'r. 


This is the general differential equation for the 
flexure of a thin plate due to uneven heating, in 
the sense that 7, 72 may be functions of both x 
and y, the variation through the plate being 
linear. P 

The functions giving the mean stress com- 
ponents X,, X, and Y, are determined by the 
equations of generalized plane thermal stress, 
and the boundary conditions, independently of. 
the flexure. 

In the problem we are considering, the tem- 
peratures 7, and 7» are functions of x only. Let 
the strip be heated and subject to the normal 
edge traction —Ea-3(7,:+T7:2). Then, by (1), 
X,=X,=0 and Y,=—Ea-3(7,+T2). The Eq. 
(5) reduces to 

1 oe 
Dy'w+Ea-—(T,+ T,)— = D-(1+¢0) —(™%-—Th2). 
2 oy" h Ox? 


- 


A particular integral is obtained by taking w to be 
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a function of x only, in which case the equation is 
satisfied provided 


Pw/ dx? =(a/h)(1+¢)(T1— 


This state of flexure corresponds to 


G,=Ih=N,=N2=0; 
Ge=(1/12)Eah(Ti1—T2). (6) 


It is maintained by bending couples (1/12)Eah? 
(T;—Te) per unit length and normal traction 
—Ea-}(7T,+T2), on the edges, together with the 
given temperature distribution. 

To make the edges free we have to superpose 
equal and opposite edge actions, with the 
stresses they produce in the interior. The effects 
of the mean normal traction Ea- }(7,+ 72) canbe 
worked out as in §3, with T7=3(7,+7.). The 
effects of edge couples —(1/12)Eah?(7T\—T>2) 
remain to be considered. 

The effects of an arbitrary distribution of 
couples on the strip of finite width could be found 
by methods analogous to those employed by 
Filon and yon Karman in the solution of the 
plane stress problem. We shall consider here only 
the case when the couples are so localized that 
one-half of the strip can be treated asa half-plane. 

The deflection of the half-plane y>0 due to a 
concentrated bending couple at the origin is** 


M 
———[2x log r—(1+0¢)y@] (7) 
~ Dx(1—o)(3-+0) 


oct 
ev 


and the deflection due to a distribution of bend- 
ing couples, G per unit length, between A and B 
(Fig. 2), obtained by integrating (7), is 


——T[r/? log-ri-sin 26; 
~ Dal ‘1 —a)(3+oe) 
— Pro" log ro° sin 262 — 2y?(01 — 2) 
+3(1 +.) (1170; — 12702) |. (8) 


The deflection due to an arbitrary distribution 
of edge couples follows from (7) or (8) by 


* A. Nadai, Elastische Platten, p. 204. 
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integration. On evaluating the couples G, and G,» 
corresponding to (8), it is found that at the edge, 
(where of course G2=G within AB, and G.=0 
outside AB), 


Gi= —(1—¢0)G/(3+c) within AB, 


and G,=0 outside AB. 
Thus, when we assemble couples on a series of 
infinitesimal segments of the edge, the couple G» 
at the edge, in any segment, is due entirely to the 
couple acting on the segment itself, and for any 
distribution of applied couple G(x) we shall have 


(G1) y-0 = G(x)(1—o) /(3+¢). 


The distribution G(x) = — (1/12) Eah?(T,—T2) 
occurring in the thermal stress problem, therefore 
gives rise to the stress couple 


G,=(1/12)Eah[1 —¢)/(3+6)](TM— 


in the edge. The extreme fiber stress correspond- 
ing to this is 


+}3Eal(1—¢)/(3+0¢) }(T1—T?). 


T») 


There is no contribution to this stress component 
from (6). The contribution (X,) from the mean 
stress is Ea-}(7,+T7>2). The extreme fiber stress 
in the edge has therefore the values 


Ea} T,+T2+[(1—0)/(3+¢) }(11—T2)}, (9) 


independent of the thickness. With localized 
heating the stress due to the temperature differ- 
ence will have the same kind of distribution in the 
interior as the mean stress, and can be visualized 
from Fig. 3, X, representing G;, and Y,, Gs. 
There will be, with 7; > 7+, a tensile stress in the 
edge, given by (9) when the positive sign is 
taken. 

As an example to illustrate heating on one face 
only, let the other face remain cool, so that T;=0. 
Then the maximum tensile stress set up is 
2EaT imax./(3+¢). For the same maximum tem- 
perature this is, for steel, about 0.61 of the stress 
EaT max. set up by heating uniformly through the 
thickness. It may thus be advantageous to keep 
one face of the strip as cool as possible. 
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Variations of the Rigidity Modulus with Changes in Applied Stress 


Guiapys M. Francis, New Jersey College for Women 
(Received December 30, 1935) 


An attempt to determine by an interferometer method whether or not any variations occurred 
in the rigidity modulus similar to those which occurred in Young’s modulus in D. K. Froman’s 
experiment. Any variations of more than 0.1 of a fringe could readily be detected. No such 


variations occurred. 





INTRODUCTION 


HE purpose of the present experiment is to 

discover whether or not any changes occur 
in the rigidity modulus when stresses are applied 
in successive steps, starting with a small stress 
applied in one direction and increasing the stress 
in the opposite direction until the elastic limit of 
the material under investigation has _ been 
reached. The elasticity is measured after each 
increase in stress by adding a small additional 
weight and observing the twist produced in the 
rod. Any changes in the rigidity modulus taking 
place would appear in the measurement of the 
twist. 

This investigation was suggested by an experi- 
ment performed by Darol K. Froman of the 
Ryerson Physical Laboratory, University of 
Chicago.!. Froman, using an_ interferometer 
method for determining the strain, showed that 
for small stresses, Hooke’s law did not hold 
when a determination of Young’s modulus was 
made, but that Young’s modulus had a maximum 
at a small but finite stress. This maximum 
occurred when the stress was approximately 
310° dynes per sq. cm. Since the constancy of 
the rigidity modulus is a proof of Hooke’s law, 
the present experiment was performed to dis- 
cover whether or not the same or similar devia- 
tions occurred when a determination of this 
modulus was made. 


EXPERIMENT 


In this experiment, the same general method 
that Froman followed in his experiment for 
determining the strain was used; namely, the 
interferometer method. 

The experimental rod L (Fig. 1) was placed 


'Froman, Phys. Rev. 35, 264 (1930). 


inside a brass tube 77. The end G of the rod L 
fitted loosely in a hole in a vertical support 
rod E, while the other end F was clamped to a 
metal piece which slid into the openings of a 
beam of a laboratory balance. The knife-edge 
of the balance, which was in a line with the 
center of the rod, rested in an agate groove which 
was clamped to the support rod Q. Both of the 
support rods were fastened rigidly to clamps on 
an optical bench so that the experimental rods 
could easily be changed without making serious 
changes in the set-up. 

At the ends of the beam of the laboratory 
balance were suspended the pans X and JY. 
The torque was applied to the rod by adding 
weights to pan X. The end G of the rod L was 
prevented from turning by means of the clamps 
I, J, K, and N joined by rods. Each time a 
different mass was placed in pan X these clamps 
were adjusted until the beam of the balance was 
level. 

Since the masses added at pan X produced a 
turning moment at N of the support £, another 
balance pan A, was added at P, the arm (ef) 
being equal to arm (cm) of the applied force. 
By adding to this pan A; at a distance (ef) 
masses equal to those on pan X a couple resulted 
so that the system balanced. 

At B, the tube 7/7 was tightly fastened to the 
rod L by a clamp which was arranged so that it 
touched the circumference of the rod in a line 
lying in a plane perpendicular to the axis of 
the rod. 

As shown in the diagram (Fig. 1) the tube // 
was not as long as the experimental rod so that 
arrangements for supporting the end C of the 
tube were necessary. This was done by fastening 
a brass bar W to the tube at the point C;. To the 
end A of the- bar a knife-edge was fastened. 
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This knife-edge rested in an agate groove in 
line with the axis of the rod. 

At C, a brass bar S was attached to the tube /7/. 
At the end U of the bar, one of the half-silvered 
mirrors VM was fastened. The other half-silvered 
mirror 4, was suspended by a bar V which was 
tightly clamped to the experimental rod L at 
the point D. This clamp was constructed so 
that it touched the rod in the same way that 


‘tthe clamp at B did. The object of these line 


contacts between the rod and clamps was to 
furnish a definite beginning and end for meas- 
uring the length of the rod. Since the arm of the 
mirror M was connected to the tube H which 
was clamped to the rod L at B, and since the 
arm of the mirror M, was attached to the rod L 
at D, any torque applied to the rod producing a 
twist between B and D tended to cause the 
mirrors to separate. 

The mirrors were placed with their centers 
10.3 cm from the center of the rod L in order to 
increase the relative change in distance between 
them when the elasticity of the rod was 
measured. ~ 

In order to prevent a torque from being 
produced by the weight of the mirror M and 
the bar S, it was necessary to arrange counter 
weights O and P;, on opposite sides of the tube 
to balance the system. The counter weight O 
was fastened at the opposite end of the bar S to 
which the mirror M was attached and was added 
to balance the tube on the knife-edges so that 
the weight of the tube was concentrated at the 
axis of the rod L. The counter weight P; was 
fastened to the bar S at the same point as the 
mirror M. This weight was added so as to bring 
the center of gravity of the mirror and bar S to 
the axis of the rod. 
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Fic. 1. Arrangement of apparatus to measure torque. 
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As a source of light a mercury are lamp, 
operated on alternating current, was used with 
filters of didymium and H. R. Yellow to cut 
out all the mercury lines with the exception of 
the green one. Since the apparatus in Fig. 1 
was very sensitive, the lamp was placed on a 
board which was suspended by springs in order 
to damp vibrations which might be transmitted 
to the apparatus. 

Between the lamp and mirrors was placed a 
large lens to convert the beam of light into a 
parallel beam. A metal sheet about a foot square 
with a circular hole of diameter one inch was 
supported between the lens and the mirrors to 
reduce the size of the beam of light so that it 
was smaller than the mirrors. 

The fringes, which were circular, were ob- 
served through a telescope, the position of which 
was kept as nearly constant as possible to ensure 
that the same point on the mirrors was observed 
each time. 

The elasticity of the rod was measured by 
applying a small additional mass to pan X and 
observing the fringe shift which it caused. The 
small mass was suspended by a ring from a hook 
above the pan of the balance. A support was 
arranged in such a way that by raising it, the 
mass would be lifted from the ring, and by 
lowering it, the mass would be added to the ring. 

It was essential that the sma!] additional mass, 


‘used for measuring the elasticity, strike the 


support quite squarely; otherwise, the mass 
would be set swinging, which in turn set the 
whole apparatus vibrating, making it impossible 
to count the number of fringes. 

In the arrangement of the present experiment 
it was possible to obtain readings for stresses 
applied in either direction since masses could be 
added to either pan. To find out what happened 
as the rod was twisted through a point when ao 
torque was applied, a small mass was added to 
pan Y and the elasticity determined for suc- 
cessive additions of masses to pan X. 

Seven rods were used in the experiment, of 
which five, the aluminum, copper, iron, steel and 
quartz rods, were solid and two, glass and steel, 
were hollow. These rods were regular commercial 
rods that are usually sold and were approxi- 
mately one meter long. The diameters averaged 
0.95 cm. 
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Fic. 2. Graph showing value of rigidity modulus as related 
to applied torque in 10° dyne cm. 


The small additional mass varied for different 
rods. For the aluminum, copper, hollow steel 
and quartz rods, it was one gram. For iron and 
solid steel rods, it was three grams, and for the 
hollow glass rod it was one-half a gram. The 
mass used depended upon the number of fringes 
which were obtained, and was such that the 
fringe count was approximately twenty-five, so 
as to make results as accurate as possible. It is 
estimated that the fringe count is accurate to 
0.1 of a fringe. 


DATA 

The distance between clamps in all cases ‘L’ 90.5 cm 
The wave-length of light used in all cases 5461 A 
Arm of applied force 15.8 cm 
Arm of mirrors (from center of tube to center of 

mirror) 10.3 cm 
Weight of ring and wire support from which sus- 

pended the small additional mass 0.681 ¢ 


DISCUSSION 


A determination of the absolute value of the 
rigidity modulus was not of primary interest in 
this experiment but it was desired to discover 
whether there were variations in the modulus 
with changes in the applied stress. It was not 
necessary in this experiment, as it would have 
been in the determination of the absolute value 
of the rigidity modulus, to have the rods uniform 
or to make exact measurements of the length of 
the rods and arm of the applied stress. These 
factors were kept constant so that if any varia- 
tions had occurred they would have appeared in 
the fringe count. The slight deviations which 
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were obtained were due to experimental errors in 
the fringe count since they appeared very 
irregularly and occurred at different points when 
measurements on the rods were repeated. 

In a determination of the rigidity modulus, 
the outer layers of the metal rod take more of 
the shearing stress than the inner layers do. 
It would be expected, therefore, that when the 
shearing stress is gradually increased the rigidity 
modulus for the outer layers would reach a 
maximum value before the inner layers would, 
provided such a maximum were to occur. This 
would mean that decreasing values of the rigidity 
modulus for the outer layers might be occurring 
at the same time that the inner layers were 
increasing to a maximum. The effect of the 
layers for which the rigidity modulus is increasing 
would thus counteract the effects of those for 
which the rigidity modulus is decreasing. There- 
fore, it would be expected that the value of the 
rigidity modulus would remain at a maximum 
value until most of the layers had been subjected 
to the torque which gives this maximum, but 
when this had occurred the value of the rigidity 
modulus would decrease to a constant value. 
However, for the solid rods this would probably 
never occur because the elastic limit of the rod 
would be reached before this would happen. 

In view of this then, if any variations were to 
occur, it might be expected that they would be 
greater for a tube than for a solid cylinder. 
It was for this reason that the tubes were used 
in this experiment. 

However, no variations in the rigidity modulus 
occurred in either the solid rods or tubes. It is 
therefore concluded that the modulus of rigidity 
is constant when the stress is gradually in- 
creased from a small stress in one direction 
through a point in the opposite direction where 
the stress is zero and then up to the elastic 
limit of the material under investigation. 

I wish to express my appreciation to the New 
Jersey College for Womea for the use of the 
laboratory and to Dr. Frank R. Pratt for his 
cooperation and helpfulness in this experiment. 
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The Synthesis of Gear Tooth Curves*: + 


J. Str_Es BeGGs, Kodak Research Laboratories, Rochester, New York 
(Received October 24, 1935) 


This paper sets forth a method of arriving analytically at a pair of gear tooth curves which 
will fulfill certain prescribed conditions. The method may also be used to find a curve which will 
mesh with a given curve and produce a constant angular velocity ratio, e.g., the mating curve 


for a circular pin tooth. 





HE analytical conditions to be met by gear 
tooth curves may be listed as follows: 


I. The primary requirement is that the two gears shall 
rotate with a constant angular velocity ratio, i.e., the pitch 
circles! shall roll on each other without slipping. 

II. Referring to Fig. 1, the point of tangency of the pitch 
circles, i.e., the pitch point, M, is the instantaneous center 
of either gear with respect to the other. 

III. The point of contact (tangency) of the tooth profiles 
at any instant is denoted by N. Since the tooth profiles are 
in rolling contact and in general, slipping, the line, MN, 
must be normal to both tooth curves at each successive 
point of contact, by condition II, otherwise, the teeth 
would either cease to make contact or else tend to pass 
through each pther (Fig. 1). 


MATHEMATICAL FORMULATION OF THE PROBLEM 


For the purpose of mathematical analysis, it 
will be convenient to think of one gear as fixed in 
the system, XOY (Fig. 2), and the other in the 
system, xoy. 

Let a+b=c 


by (1) a#=b¢; therefore, (1) 


(2) 


The equations of transformation between the 
two systems are as follows: 


6+o=c/b-0. 


x=(Y—c sin 6) sin c/b0 


+(X—c cos @) cosc/b@, (3) 
y=(Y—c sin @) cos c/bé 
—(X—c cos @) sinc/bé. (4) 


Let the equation of the tooth curve of the gear 
in XOY (Fig. 1) be given by 


Y=F(X). (5) 





*A revision of a paper read before the Mathematics 
Club of the University of Alberta. 

+ Communication No. 560 from the Kodak Research 
Laboratories. 

‘For the meaning of technical terms, see, e.g., Machin- 
ery’s Handbook (Industrial Press, New York, 1934). 
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Then, by condition III, the point of contact, 
N, lies at the foot of the perpendicular from M/ 
to the curve of Eq. (5). 

Obviously, N will trace out in xoy the pinion 
tooth curve. Let this be given by 


y=f(x). (6) 


Eq. (6) may be found by taking the co- 
ordinates (X, Y), of N and transforming them by 
(3) and (4) to the system xoy, as paramétric 
equations in 0. 

The equation of any straight line through / 
is given by 

Y=a sin 0+k(X —a cos @). (7) 


Since it is perpendicular to Y= F(X), at N, by 





\ vere) 


Fic. 1. The general form of gear tooth curve. M is the 
pitch point and N the point of contact of the curves. 
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Fic. 2. The gear is considered as fixed in the coordinate 
system XOY and the pinion in the system xoy. 
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Fic. 3. The form of tooth which gives a constant bearing 
load at constant torque. This is found to be the 
involute. 


condition ITT, 
k=—1,dF(X)/dX, at the point N. (8) 
Therefore, (7) becomes 
Y=a sin 6—1/ F’(X)(X —a cos 6). (9) 


The coordinates of N are given by the simul- 
taneous solution of (5) and (9). 


THE INVOLUTE TooTH 


As an example of this method, we impose the 
following conditions: The reaction on the bear- 
ings is to be constant for a constant torque, since 
a variation in this force would tend to set up 
vibration and noise. This means that the normal 
JN must make a constant angle, a, for instance, 
with the line of centers of the gears. 

Let Eq. (5) of the required gear tooth curve be 
given by 





X¥ = X(6), (10) 
Y= Y(@). (11) 
From Fig. 3, 
dY dX = tan (@+a). (12) 
Therefore, (7) becomes 
X —a cos 6 
Y=a sin @—— --, (13) 
tan (6@+a 
Then, differentiating (13), 
dY=a cos 6d6+ (X —a cos @) 
x sec? (0+a) cot? (@+a)d0— cot (@+a)dX 
—asin 6 cot (@+a)d@. (14) 


The simultaneous solution of (12) and (14) 
gives the differential equation of the locus of N. 
Substituting from (12) in (14) for dY, we have 
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Fic. 4. The involute curve. 
dX cos (@+a) cos? (@+a) 
Hn COG Ea—————————, (55) 
dé sin (0+ a) sin (@+a) 


which is a linear differential equation whose 
solution is 


X =a cos af (6+ tan a) sin (@+a) 


+ cos (@+a)]. (16)? 
Similarly, Y=acos a[sin (@+a) 
—(6+ tan a) cos (@+a) ]. (17) 


Referring to Fig. 4, the general equation of the 
involute is given by 


X =d cos ¥+d(y¥+8) sin y, (18) 

Y=d sin Y—d(W+8) cos y. (19) 
Comparing (16) with (18), let 

a cos a=d, (20) 

6+a=y, (21) 

a+B= tan a=(a?—d?)!/d. (22) 


The curve (16) and (17) is, therefore, an invo- 
lute with the base circle, a cos a, and starting at 
an angle, a= tana—a, below the X axis, the 
constants being so chosen as to make the pitch 
point lie on the X axis. 

By imposing certain initial conditions, we have 
thus arrived at a unique solution of our problem, 
i.e., the involute is the only curve which will 
satisfy these conditions. 

To find the mating tooth curve, we transform 
Eqs. (18) and (19) by means of (3) and (4), 
obtaining, 


x=-—b a-d[cos (—a/b@+a)+(—a,/bé 


+tan a) sin (—a/b@+a)]}, (23) 





2T am indebted to Miss Helen Shea of the University 
of Toronto for the formal integration of Eq. (15). 
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y= —b/a-d[sin (—a/b0+a) —(—a/b0 


+tan a) cos (—a/b0+a)}. (24) 


The mating curve is also seen to be an involute. 

As slipping is an important factor in the wear 
of gear teeth, it is desirable to have an expression 
for it. Referring to Fig. 1, the slip is seen to be 
equal to JN multiplied by the relative angular 
movement and is, therefore, zero, only when the 
point of contact coincides with the pitch point. 

Let dS be an element of arc of the involute gear 
tooth, and ds be an element of arc of the pinion 
tooth. If o represents the slip, then 


do=dS—ds=d(1+a/b)6d0. (25) 
do /d0@=d(1+a/b)é. 


id d-c 
p=otoo= dao=—#*. (27) 
6=0 2b 


The rate of slipping (26) 


Slip from 


In practice, to facilitate the cutting of the 
teeth, standards have been set up. The pressure 
angle, a, is taken as 143° and 20°. From Eqs. (18) 
and (19), we see that the angle of obliquity and 
pitch diameter determine the shape of the tooth. 
The diametral pitch determines only its size. 


MATING CURVES 


This method may also be used to find the 
mating curve to any given curve, (5), e.g., 
suppose we wish to find the shape of tooth which 
will mate with a circular pin tooth and satisfy 
condition I. In this case, (5) is of the form 
(Fig. 5): 
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Fic. 5. The pin tooth and its mating curve. 


(28) 


(29) 


X =a-—d sin 6/2, 
Y=d cos 6/2. 


Transforming by Eqs. (3) and (4), we obtain 
the mating curve 


x=acosc/b@+d sin (1/2+a/b)0 


—ccosa/b@ (30) 
y=—asin c/b0+d cos (1/2+a/b)0 
+csina/b@, (31) 


which is not a circle. 

It is interesting to note that if we replace the 
pin by a ball bearing of the same diameter and 
thus introduce at least three balls and an outer 
race unfixed in either YOY or xoy, we obtain a 
pair of gears in which there is no sliding friction. 

Any form of (5), such as an involute or cycloid, 
which will undergo the transformation (3) and (4) 
without change of form will yield a mating tooth 
of similar shape, or, if the gear and pinion are of 
equal diameter, of identical shape, which is an 
obvious advantage from the manufacturing 
standpoint. 
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LETTERS TO THE EDITOR 


This section will accept brief reports of new work 
of general interest to industrial physicists. The 
Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Con- 


A Formula for the Birefringence of Vibrating Media 


Consider an optical arrangement (such as used by Giin- 
ther') in which a transparent medium is placed between a 
pair of crossed Nicols and is traversed by a beam of parallel 
monochromatic light. The diagonals of the Nicols are set 
at 45° with respect to the principal plane of the medium. 
When the medium is caused to vibrate, the field of vision 
brightens and/or darkens in those regions in which the 
mechanical strain is sufficiently great. The following theory, 
which has been used by the writers for the past two years, 
permits a quantitative analysis of the vibrations of trans- 
parent media, as for example piezoelectric crystals. 

The intensity J of the light transmitted through both 
Nicols is? 

T=TI max sin? tR, (1) 
where R) is the retardation of the extraordinary ray with 
respect to the ordinary ray. If the medium be subjected 


to a force field F sin (27t/T), one can show from Pockels’ 
theories* for the linear piezo- and electro-optical effects that 


Ry =cLF sin (2xt/T), (2) 


where L is the thickness of the medium and ¢ is a linear 
function of the optical constants. Combining Eqs. (1) and 
(2) one now finds the instantaneous intensity to be 


1; =Imax sin? (xcLF sin (22t/T)) 
= (Imax/2)[1—cos (24cLF sin 22t/T) ] 
= (Imax/2)[1—cos (x sin @) ]. (3) 


When the frequency 1/7 is more rapid than 20 cycles per 
second, the eye discerns the average intensity J,. 


Imax 2 rh . 
l= fy [1—cos (x sin @) ]dt 
Ima: °r . ‘ 
a | [1—cos (x sin @) ]d@ 
2r v0 


tributions to this section must reach the office of 
the Managing Editor not later than the twenty-fifth 
of the month preceding that of the issue in which the 
letter is to appear. 


Tmax Tinax 
= [1—Jo(x) }#=— eh — Jy(2ncLF) 1. (4) 


As F, and therefore x, is increased, J, increases until it 
reaches a maximum at x =3.83. With increasing values of 
x, I, decreases to a minimum at x=7.02, reaches another 
maximum at x=10.17, a second minimum at x=13.32, 
etc. For every observed maximum or minimum the corre- 
sponding value of x is thus known. 

This theory explains the beautiful photographs by 
Briininghaus.® It also explains his ‘‘polarization chroma- 
tique”” when one includes the variation of ¢ with wave- 
length. 

The first minimum of J, Eq. (1), is reached in the static 
case when R,=1 or F=1/cL. The corresponding minimum 
is reached in the dynamic case when F=(7.02/27)(1/cL) 
=1.12/cL, or at a value 12 percent greater than in the 
static case. 

For values of x less than 1, J, is given with an accuracy 
greater than 93 percent by 


Ia =Imaxx?/8. 


This equation predicts the intensity of light which passes 
through the Okolicsanyi cell,®»7 an instrument used in 
television. 

Papers relating to the application of this method to 
special problems in photoelasticity will be published else- 
where. 

Joun W. Cookson 


HAROLD OSTERBERG 
University of Wisconsin, 
Madison, Wisconsin, 
March 25, 1936. 
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